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Abstract. 

The implications of the WMAP (Wilkinson Microwave Anisotropy Probe) third 
year data for inflation arc investigated using both the slow-roll approximation 
and an exact numerical integration of the inflationary power spectra including a 
phenomenological modelling of the reheating era. At slow-roll leading order, the 
constraints t\ < 0.022 and —0.07 < €2 < 0.07 are obtained at 95% CL (Confidence 
Level) implying a tensor-to-scalar ratio r\o < 0.21 and a Hubble parameter during 
inflation H/m Pl < 1.3 x 10~ 5 . At next-to-leading order, a tendency for £3 > is 
observed. With regards to the exact numerical integration, large field models, V{<j>) oc 
cj) p , with p > 3.1 are now excluded at 95% CL. Small field models, V(<p) oc 1 — (cf>/ fi) p , 
are still compatible with the data for all values of p. However, if /i/m P1 < 10 is 
assumed, then the case p = 2 is slightly disfavoured. In addition, mild constraints on 
the reheating temperature for an extreme equation of state u> ro h ^ — 1/3 are found, 
namely T rc h > 2 TeV at 95% CL. Hybrid models are disfavoured by the data, the 
best fit model having A\ 2 — +5 with two extra parameters in comparison with large 
field models. Running mass models remain compatible, but no prior independent 
constraints can be obtained. Finally, superimposed oscillations of trans-Planckian 
origin are studied. The vanilla slow-roll model is still the most probable one. However, 
the overall statistical weight in favour of superimposed oscillations has increased in 
comparison with the WMAP first year data, the amplitude of the oscillations satisfying 
2 1 a; I tr < 0.76 at 95% CL. The best fit model leads to an improvement of A% 2 ~ — 12 for 
3 extra parameters. Moreover, compared to other oscillatory patterns, the logarithmic 
shape is favoured. 



PACS numbers: 98.80.Cq, 98.70.Vc 
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1. Introduction 

The recent release of the three years WMAP data El El E] constitutes an important 
step for the theory of inflation. One now has at our disposal high accuracy data that 
can be used to probe the details of the inflationary scenario and to learn about the 
physical conditions that prevailed in the very early universe, at very high energies 
comparable to the Grand Unified Theory (GUT) scale. There are many aspects that 
would be interesting to study but, clearly, in a first approach, one can restrict ourselves 
to (effective) single field models and see whether it is already possible to constrain the 
shape of the inflaton potential V{4>). In particular, it is important to known whether it 
is necessary to go beyond a simple Harrison-Zeldovitch (scale-invariant) power spectrum 
to correctly fit the data. If so, as indicated by the results of reference [2] and from a 
model building point of view, this means that the inflaton potential is not completely flat 
and/or that inflation is not driven by a pure cosmological constant. Equivalently, this 
also means that the observations start feeling the non-trivial shape of the potential. In 
this case, a non-vanishing second order derivative of the potential is seen while the first 
order derivative of V(<p) presumably remains unconstrained from below since, otherwise, 
primordial gravitational waves would have been detected. Besides the previous issue, one 
would also like to go further and to study which inflationary models remain compatible 
with the data and which ones are ruled out. Addressing these questions is the main goal 
of this article. Notice that we restrict our considerations to WMAP3 data only although 
including other data sets could allow us to obtain tighter constraints on the inflationary 
models studied here. In a first step, we think it is more reasonable to proceed this way 
in order not to mix the effects of using different data sets with those originating from 
the new numerical approach introduced in this article. 

To deal with this problem we proceed as follow. We first use the slow-roll 
approximation and derive the constraints put by the third year WMAP data on the 
first three slow-roll parameters (i.e. second order slow- roll approximation). Then, we 
compute exactly these parameters for the large field, small field, hybrid and running 
mass models. In particular, we show that some widely used approximate expressions 
for e\ and are no longer sufficient to assess the likelihood of some models (small 
field models) given the quality of the data. For this reason, in this article, we always 
evaluate exactly the slow-roll parameters, using simple numerical methods if necessary. 
The implications of the WMAP3 data using a slow-roll prior have also been investigated 
in references EH 13 El E] • 

In a second step, one frees ourselves from any approximation (except the linear 
theory of cosmological perturbations) and calculate the power spectra exactly by means 
of numerical computations. The models that we study are the same than the ones 
already mentioned for the slow-roll case, namely large field, small field, hybrid and 
running mass models. These exact power spectra are computed mode by mode and 
fed into a modified Cosmic Background Microwave (CMB) code, here a modified CAMB 
code [TOl, which allows us to determine the temperature and polarisation multipole 
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moments. Finally, we explore the corresponding parameters space by using Monte- 
Carlo techniques as implemented in the COSMOMC code JT] together with the likelihood 
code developed by the WMAP team [2]. This allows us to put constraints on the free 
parameters characterising the models during the inflationary phase but also, in principle, 
during the reheating phase although, most of the time, the accuracy of the data is not 
sufficient to obtain relevant limits on the reheating temperature. 

This article is organised as follows. In section |21 we briefly recall some basic facts 
about inflation, reheating and the theory of inflationary cosmological perturbations 
of quantum-mechanical origin. Section |3] is devoted to the slow-roll approximation. 
In subsection 13 .1| the derivation of the slow-roll scalar and tensor power spectra is 
recalled. In subsection 13.21 we present the WMAP data constraints on the slow-roll 
parameters at first order (i.e. e\ and €2) and also at second order (i.e. the two previous 
ones plus £3). In subsections 13 . 31 to 13 . fi[ we calculate the slow-roll parameters for the large 
field, small field, hybrid and running mass inflationary models. These results are then 
compared to the constraints on ei, €2 and €3 obtained previously. Section |U is devoted 
to the the exact computations of the inflationary power spectra for the four models 
mentioned above. In subsection 14. 1| we briefly describe the method and the code used 
to perform the numerical calculations. In subsections 14.21 to 14.51 the exact numerical 
results are used to discuss the constraints put by the third year WMAP data on the free 
parameters describing the models but also (when possible) on the subsequent reheating 
phase. Finally, in section we investigate the presence of superimposed oscillations 
in the CMB multipoles. In subsection 15. 1[ we discuss a possible physical origin for 
those oscillations, namely trans-Planckian effects during inflation. This allows us to 
use a well-motivated and well-defined shape for the oscillatory power spectra. Then, in 
subsection 15.21 we compare these spectra to the third year WMAP data and use them 
to put constraints on the amplitude, the frequency and the phase of the superimposed 
oscillations. In the last section El we recap our findings and present our conclusions. 

2. Basics of inflation 

2.1. The background 

2.1.1. The accelerated phase In this section, in order to describe the general setting 
and to fix our notations, we recall some basic and well-known facts about inflation 
at the background and perturbed levels. A phase of inflation is supposed to make the 
universe homogeneous and isotropic on large scales. It also drastically reduces its spatial 
curvature in agreement with the observations which indicate that the spatial sections 
are extremely flat. As a consequence, the metric tensor which describes the geometry of 
the Universe can be taken of the Friedman-Lemaitre-Robertson- Walker (FLRW) form, 
namely 

ds 2 = -dt 2 + a 2 (t)5 ij dx i dx j = a 2 (r]) (-dr/ 2 + <%dxW) , (1) 
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where 5ij is the Kronecker symbol. The variable t is the cosmic time while r\ is the 
conformal time and these two quantities are related by dt = adrj. In the following, both 
will be used, the choice of using one rather than the other being made for convenience 
only and according to the problem at hand. Another important time variable is the 
number of e-folds defined by the following expression 

where a- m is the value of the scale factor at some initial time. As required to solve the 
flatness and homogeneity issue of the FLRW model, the total number of e-folds during 
inflation must be greater than 60. 

The evolution of a(r]), the only free function in the above metric element, is 
controlled by the Einstein equations. If matter is assumed to be a perfect fluid, they 
read 

2 
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where a prime denotes a derivative with respect to conformal time and p, P are 
respectively the energy density and the pressure of the cosmological fluid driving the 
dynamics of the universe. In the following, we will use the conformal Hubble parameter 
defined by 7i = a' /a. In addition, one also has the energy conservation equation 

p + m{p + P) = , (4) 

which can also be obtained from the Einstein equations by means of the Bianchi 
identities. 

By definition, inflation is a phase of accelerated expansion for which the scale factor 
satisfies O Q3J d 113 

^>0, (5) 
and this condition can also be re-written as 

£lW , 1 _(^(0<, (6) 

The quantity €\ is in fact nothing but the first Hubble-flow (or slow-roll) parameter 
—H/H 2 where H = Ti/a denotes the physical Hubble parameter (see below). Inflation 
stops when e\ = 1. Another way to express the acceleration of the scale factor is to 
combine equations ©• One gets 

~ = -^(p + 3P), (7) 

where a dot denotes a derivative with respect to the cosmic time t. Clearly, inflation 
can be obtained if the fluid dominating the universe has a negative pressure such that 
P < -p/3. 

A possible implementation of the inflationary scenario is to assume that the matter 
content of the universe is described by a scalar field <p{rj) with a potential V(4>) 
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since, when V{<f>) is sufficiently flat, the effective pressure of the scalar field can be 
negative O E3 HEl UH UH] • In this article, we restrict ourselves to the case of a single 
scalar field. In this case, the two Einstein equations take the form 
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Using the expressions of p and P for a scalar field, the conservation equation (J3J) reduces 
to the Klein-Gordon equation written in a FLRW background, namely 



b" + 2H<p' + a 



2 dV(0) 







(10) 



This equation can also be directly derived from (jHJ) and 0. In fact, the physical 
interpretation is made easier if the Klein-Gordon equation is written in terms of the 
number of e-folds, namely ^7] 



d 2 (f) d(p m^dln^O) 



+ 



3-eidiV 2 ' dN 8tt 
We see that studying the evolution of the scalar field is in fact equivalent to study the 
motion of particle in an effective potential In V(<f>) with a slightly variable mass and 
a constant friction term. This is the reason why, in the following, when we consider 
concrete models, we will also pay attention to the logarithm of the potential in order to 
gain intuition about how inflation proceeds. 

If we are given a model, that is to say a concrete form for the potential V(4>), 
then the equations of motion can be integrated from the initial conditions <ft- m and 
(p[ n . In general, the shape of the potentials that we consider does not allow simple 
analytical solutions of the Einstein equations. In that case, one has to rely on analytical 
approximations of numerical calculations. However, the following potential ^H] 



V{(f>) = M 4 exp 



(12) 



is an example where the exact integration can be performed analytically. One obtains 
power-law inflation for which the scale factor and the scalar field are respectively given 
by 



a(rj) = £ \r) 



1+/9 



+ 



2 + 13 



[1 + 13) ln\r)\ 



(13) 



n(l + (3) 

In this model the parameter e\ is given by e x = (2 + /3)/(l + (3) and, therefore, inflation 
occurs if /3 < —2 (we do not consider the case where —2 < (3 < —1 which cannot be 
realised with a single scalar field). However, this model is not a satisfactory model, since 
inflation never stops, at least if one does not use another mechanism. 
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2.1.2. The end of inflation When an exact integration is performed, it is very 
important to have a description of the reheating phase. Without such a description, 
one cannot relate the physical scales today to the physical scales during inflation simply 
because one does not know the entire history of the universe. As mentioned before, 
inflation ends for <fi = en d, or p = p cn d (or again for H = iJ en d), when e\ — 1. The 
difficulty is that the process of reheating can be model dependent. For the large field 
inflation models, one may use a phenomenological description based on the fact that the 
potential is given by V(4>) oc <p p . In this case, after the end of inflation, the field starts 
oscillating around its minimum [TTH I2U|. In this regime, one can show that the average 
energy density behaves as p m { oc a~ 6p /( p+2 ) and, as a result, the scale factor is given by 
a(t) oc t(p+ 2 )/( 3 P). In other words, the equation of state P = wp during reheating has a 
constant state parameter given by 

p - 2 

^ e h = — • (14) 

For a massive scalar field, p = 2, the energy density evolves as in a matter-dominated 
epoch while for a quartic potential, p — 4, the energy density behaves as in radiation- 
dominated era. Nevertheless, one always has 

1 . . 

- - < w rch < 1 , (15) 

for p > 1 . 

The previous considerations and the description of the oscillatory phase were just 
based on the Klein-Gordon equation (|10j). It is clear that, with this equation only, the 
decay of the inflaton field into radiation cannot be described. In order to take into 
account this effect, it is common to phenomenologically add a friction term T<f> in the 
Klein-Gordon equation which now reads |19| 

+ 3#0 + T0+ ^yP~ = 0. (16) 



In fact, this description cannot account for the complexity of the reheating stage [21 
I2H 1221 123 1221231 H3 123 I2E1 121 EDI- But for our purpose, the above treatment will be 
sufficient. Since we have modified the Klein-Gordon equation one should also modify 
the conservation equation for the radiation energy density in order to ensure the total 
energy conservation: 

Prad = -4-ffp ra d + Fp inf . (17) 

In fact, the previous discussion can be generalized to other inflationary models by 
simply assuming that the reheating era takes place with a constant state parameter 
w re h for which a oc £ 2 /( 3+3w roh) anc [ p ^ a ~ 3 ~ 3Wrch . Note however that unlike for the large 
field models, this parameter is not necessarily related to the inflationary potential and 
equation ([Hjl is generally not satisfied. However, the equations for p in f and p ra d can be 
exactly integrated and the solutions read ^1] 

/ \ -3(l+w rc h) 

Pini(t) = Pend ( exp [- (1 + W vch ) T(t- tend)] , (18) 

V a end / 
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Figure 1. The quantity p r ad/Pcnd for w rc h = 2 or p = 2, t cn d = 1 and T = 0.1. The 
solid black line is the exact expression obtained from equation Ijl9|l while the dotted 
red line is the approximation of H20(l . At t/t cn< \ ^> 1, the difference between the two 
curves is approximately a factor two. Here, the time t and T are measured in units of 
m Pl . 
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.3 + 3w reh '^ ' — — 1 f ■ (19) 

where t = tend is the time at which the oscillations start (i.e. the time at which inflation 
ends). The function j(a,x) is the incomplete gamma function jHU E2j defined by 
j(a, x) = J* e~*t a_1 dt. It was shown in Ref. jEj that the expression of p ra d can be 
simplified if one uses the expression of the incomplete gamma function for small values 
of its argument. For times t > t en d with t end <C t re h where t rc h = r _1 , one obtains the 
solution [19J 



Prad(t) - np end 



^end \ 3 + 3u> re h 
t ) 5 - 3Wreh 



tend 



- (5 -3u) re h ) / (3+3«) re h ) 



(20) 



Let us notice that the quantity (5 — 3w rc h)/(3 + 3w rc h) is always positive. Therefore, 
for t 3> tend, the second term in the above equation becomes negligible. The exact (fT^j) 
and approximate (|2T?|) expressions of p m d(£) are presented in figured] At t = t re h, one 
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can consider that the phase of reheating has been completed. After thermalisation, the 
energy density of radiation takes the form p rad = g*7r 2 T 4 / 30 where g* is the number of 
relativistic degrees of freedom. Expressed at t — t re h, this quantity must be equal to 
p ra d in equation (|2Uj) and one arrives at 

1/4 

y/T^, (21) 
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where we have used p en( j = 3H 2 ad m 2 J(8ir) and t end = 2#- 1 d /[3(l + w Ieh )}. Moreover, 
The number of e-foldings during reheating can be evaluated as 

iV reh ^ In ( 3+^*^ . (22) 

3 + 3w reh V 2 r J K J 

This result is important because it allows us to relate the physical scales now to the 
physical scales during inflation. To do so, one must take into account the fact that 
a large number of e-folds can exist during reheating. This also means that the link 
between what happens during inflation to what happens now depends on the details 
of the reheating period. In our simplified model, this period is described by a single 
parameter T or, equivalently, N^. One can also combine the two above equations in 
order to express the reheating temperature only in terms of N reh and u> re h- This leads 
to 

1/4 

e -3(l+«; rc h)A r r oh/4 ^3) 



1/4 30 1/4 1/4 f 3 + 3w ieh \ 



This expression also depends on p end but this quantity is known once the background 
inflationary evolution is solved. 



2.2. The co smological perturbations 

Having integrated the background, we now turn to the cosmological perturbations 
themselves. The power spectrum can be obtained by a mode by mode integration. 
Before addressing this issue, let us recall some basic definitions. The perturbed line 
element around a spatially flat FLRW universe can be expressed as [3*3*| I34j : 

ds 2 = a 2 (t]) {- (1 - 2(f)) dry 2 + 2d i B&x i &r) + [(1 - 2ip) 5y 

+ 2d i d j E + h lj \<\x i dx j } . (24) 

The four functions <ft, B, if) and E represent the scalar sector whereas the transverse 
and traceless tensor hij, i.e. satisfying hi = $hij = 0, represents the gravitational 
waves. There are no vector perturbations because a single scalar field cannot seed curly 
perturbations. At the linear level, the two types of perturbations decouple and therefore 
can be treated separately. Since the scalar sector suffers from a gauge dependence, it is 
more convenient to work with the gauge-invariant Bardeen potentials $ and the gauge- 
invariant perturbed scalar field S(f)^ defined by |HSl E3j 

$ = + - [ a (B - E'))' , 5<f) igi) = 5(f) + 4>' (B — E') . (25) 
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In the first of the two above equations, <p denotes the perturbed time-time component 
of the metric element (J2*%jl an d not the perturbed scalar field. In fact, it is clear that 
one has only one degree of freedom because $ and S(p^ are related by the perturbed 
Einstein equations. One can therefore reduce the study of the scalar sector to the study 
of a single variable, namely the so-called Mukhanov- Sasaki variable defined by 



!,, +4 



(26) 



In order to set quantum initial conditions, it turns out to be more convenient to work 
with the rescaled variable fi s = —y/2K,v. On the other hand, the quantity of interest in 
the primordial power spectrum is the comoving curvature perturbation ( (rj, x) defined 
by 

C = -^5=- (27) 



Since no tensorial quantity can be used to generate an infinitesimal coordinate 
transformation, the tensor sector is gauge invariant and fully characterised by the 
quantity /i T related to through the relation 

hij = —Qij ■ (28) 
a 

The are the transverse and traceless eigentensors of the Laplace operator on space- 
like hypersurfaces. 

The central result of the theory of inflationary cosmological perturbations is that 
the quantities /z SjT (or rather their corresponding Fourier transform) both obey the 
equation of motion of a parametric oscillator [HEl ESI E3 EB] 

^ + Ws, t (MKt = 0, (29) 

where the time variation of the frequencies only depend on the behaviour of the 
background and is given by 



ul(k,r 1 ) = k 2 - uUk,r}) = k 2 - — . (30) 

aJe{ a 



In these expressions k denotes the comoving wavenumber. In order to solve the above 
equations, one postulates that the quantum fields are initially placed in the vacuum 
state when the mode k is well within the Hubble radius, which amounts to assume that 

lim tisT(V) = T 1= — i (31) 

k/(aH)^+oo^ ' TOpi J2k 



where 7/; is the initial conformal time at the beginning of inflation. Then, a mode by 
mode integration of equation ()29|) allows the determination of /z SjT for a given k at any 
time rj. 

It is then straightforward to determine the resulting power spectra. From a 
calculation of the two-point correlation functions, one obtains 
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These are the quantities that seed the subsequent CMB anisotropies. Usually, the 
properties of these primordial power spectra are characterised by the spectral indices 
and their "running" . They are defined by the coefficients of Taylor expansions of the 
power spectra with respect to Ink, evaluated at an arbitrary pivot scale fc*, namely 

dln{k 3 P h ) 



dln( k 3 P c ) 

71= - 1 = 



din A; 



fin 



d In k 



(33) 



For the runnings, one similarly has the two following expressions 
_ d 2 ln(PP c ) _ d 2 \n(k 3 P h ) 

a -^J^ t ; aT =^o^r t ,- (34) 

and in principle, we could also define the running of the running and so on. 

3. Testing the slow-roll models 

3.1. Basics of slow-roll inflation 

It was shown in the previous section that the knowledge of the background evolution 
is sufficient to calculate the time-dependent frequencies u; s (k, rj) and u> T (k,rj) which are 
the only quantities needed to integrate (possibly numerically) the equations of motion 
of the cosmological perturbations. The slow-roll approximation allows to perturbatively 
estimate ui BjT (k, rj) and, hence, to derive approximated expression for the scalar and 
tensor power spectra (231101 E]- Although there are several definitions of the slow-roll 
parameters, in this article, we choose to work with the Hubble-flow parameters {e„} 
defined by 



d In I e n I 

en+i = dN , n > , (35) 

where, as already mentioned, N is the number of e-folds since some initial time r]i n . The 
above hierarchy starts from eo = H m /H. With this definition, all the e n are typically 
of the same order of magnitude. One has slow-roll inflation as long as \e n \ <C 1, for all 
n > while, as already mentioned, inflation takes place if ex < 1. 

The physical interpretation of these parameters has been discussed in Ref. Let 
us briefly recall the results for e\ and e<i- Although the e n parameters make no reference 
to the matter content of the Universe (they are only defined in terms of the expansion 
rate), it is nevertheless interesting, when we assume that a scalar field <fi is responsible 
for inflation, to express them in terms of cf). One obtains [13] 



t,= V/2 + V(*)' d;(Yj = M2- f ')- (36) 

From the above expressions, one sees that ei/3 measures the ratio of the kinetic energy 
to the total energy (i.e. kinetic plus potential energy). Using the link between e\ and 
e 2 , one has 

ei = H ei e 2 . (37) 
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Therefore, given the fact that e\ is positive definite, e 2 > (respectively e 2 < 0) 
represents a model where the kinetic energy itself increases (respectively decreases) with 
respect to the total energy. From equation (JHEJ), ^2 = 2ei marks the frontier between 
models where the kinetic energy increases (e 2 > 2ei) and the models where it decreases 
(e 2 < 2ei). 

It was demonstrated in |4*3*1 H3] that the Hubble-flow parameters can be expressed 
in terms of the inflaton potential and its derivatives. The exact expressions read 

m = *S*<tl-<±Y (38) 



8tt V 3 



- m ^tv-i+ih (39) 

^ = 3tf 2 f 2 Cl - ^ - -c? + -e l£2 - -e 2 - le*^ . (40) 
d0 2 V 2 3 6 12 6 J K ' 

In order to have a meaningful model, the potential should be positive and, hence, e\ < 3. 

In fact, it is more interesting to express the slow- roll parameters in terms of the potential 

and its derivatives because, in practice, a model is defined by its potential. Therefore, 

one has to invert the above expressions. At leading order in these parameters, one 

obtains 



8tt _ m 2 {V'\ m 2 
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(42) 



where, in the present context, a prime denotes a derivative with respect to the scalar 

field (f>. Since we have pushed the calculation up to order two in the slow-roll parameters, 

one may be worried about the fact that we only perform the inversion at leading order. 

Following [31)] > we can define two new slow- roll parameters 

,2 /y>\ 2 m 2 yli 

Vv = ^T7- (43) 



16tt \ V J ' 8tt V 

At the next-to-leading order, one obtains 

4 2 

ei = ev - g4 + 3 £ v^v + • • • • (44) 

At leading order, e\ = e v and one recovers the expression (}4~Tj) for the first Hubble-flow 
parameter. As expected, the corrections to this expression are quadratic in the slow-roll 
parameters. Therefore, as long as the slow-roll parameters are small, the predictions 
for a concrete model, i.e. the corresponding location in the (ei,e 2 ) plane, are almost 
unchanged and, in any case, undetectable with the current data. We will come back 
to this issue in the following and check this claim explicitly for large fields models. Of 
course, if one wishes to determine, say, the scalar spectral index at quadratic order, that 
is to say to include the quadratic terms in the expression of n s — 1 (see below), then it 
would be mandatory to perform the inversion at quadratic order as well. 



Inflation after WMAP3 



12 



We now turn to the perturbative expression of the power spectra. The strategy 
consists in expanding the power spectra about the pivot scale fc*. The choice of this 
particular wavenumber is a priori arbitrary but must be chosen in a way minimising 
the uncertainties coming from the perturbative expansion. Therefore, a good choice is 
around the middle of the range of scales probed by the CMB (the size of which is about 
three decades). Consequently, the usual choice is k* = 0.05 Mpc -1 . The expression of 
P{k) reads 

^ - ao + ai Jf] + ^(h\ + ..., (45, 



kPPo^k*) \k* J 2 V k 



where 



2 



k 3 P Co (h) = , k 3 P ho (K) = 16-^- . (46) 

7reimpj nm^ 

The coefficients a« are then determined in terms of the slow-roll parameters. These 
calculations are non trivial especially when one goes beyond the leading order because, 
in that case, the usual approximations in terms of Bessel functions are no longer valid. 
It is therefore necessary to use more sophisticated methods either based on the WKB 
approximation |4^ Wf\ |4~9~] and/or on the Green function methods [SUJ EIJ. The 
latter method has been used in several works (5UJ Elj to derive the coefficients a 4 and 
we only quote the results. For scalar perturbations, one gets 

f 2 

2 



af = 1 - 2 (C + 1) e a - Ce 2 + ( 2C 2 + 2C + — - 5 ) e\ 



> 2 + Qe 2e 3, (47) 

af = - 2ei - e 2 + 2(2C + l)e\ + (2C - \)e x e 2 + Ct\ - Ce 2 e 3 , (48) 
4 S) = 4e 2 + 2e ie2 + e 2 - e 2 e 3 , (49) 

where C = 7 E + ln2 — 2^ —0.7296, 7 E being the Euler constant. For the gravitational 
waves, the coefficients a« read 

(T) _ 1 o (n , -, \ r , or y2 , o^f , " c L2 



= 1 - 2 (C + 1) ei + \2C l + 2C + — - 5 j ef 

+ - 2C + ^ - 2^ Cl e 2 , (50) 

a« T) = - 2ei + 2(2C + l)e 2 - 2(C + l)e x e 2 , (51) 
4 T) = 4e 2 - 2eie 2 . (52) 

We see that each coefficients a« starts at order e l n . In fact, this can be proven analytically 
since the coefficients a^i can be obtained from the coefficient a« by differentiation with 
respect to the number of e-folds. Equivalently, if one is able to determine the coefficient 
a at order e % n , then one can derive a non- vanishing expression for the coefficients a,j up to 
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j = i. The above considerations also lead to the domain of validity of the expression (J45j) . 
Clearly, one has a meaningful expansion as long as 

e n ln^J<l, (53) 

where we have again used the fact that is of order e l n . In principle, this estimate is 
valid for any slow-roll parameters. Given the fact that the CMB observations probe 
about three decades in wavenumbers and assuming that a natural location for the pivot 
scale has indeed been chosen in the middle of this range, one arrives at 

|e„|<0.29. (54) 

In practice, in the rest of this article, we will consider that the equation (|4~5*|) remains 
valid as long as |e n | < 0.1. 

The logarithm of the power spectra can also be Taylor expanded in terms of the 
logarithm of the comoving wave number: 

n k 3 P(k) , , , fk\ &2, o / k\ 

where the coefficients hi are given by 

6« S) = - 2 (C + 1) e x - Ct 2 + (-2C + y - 7^ e{ 

+ (-^-3C + ^-7) ei e 2+ (^-l)^ 

+ f4c 2 + ^V 2 e3, (56) 



2 24, 

hf = - 2ei - e 2 - 2el - (2C + 3) ei e 2 - Ce 2 e 3 , (57) 
b w = _ 2eie2 -e 2 e 3 , (58) 

for the scalar power spectrum, while the b^s associated with the tensor power spectrum 
read 

b^= _ 2 (C+l)e 1 + (-2C + ^-7 



2 



-C 2 -2C+^-2j ei e 2 , (59) 

C = - 2ei - 2e? - 2(C + l) ei e 2 , (60) 

6 2 T) = - 2eie 2 . (61) 

The interest of the coefficients bi is that they are directly related to observable quantities, 
namely to the spectral indices and the runnings. Explicitly, one has 

n s -l = bf\ n T = b ( P, a s = b 2 s) , a T = b 2 T) . (62) 

In particular, one sees that, at first order in a consistent slow-roll expansion, the scalar 
and tensor runnings vanish. 



2 
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Finally, the ratio of amplitudes of scalars and tensors at the pivot reads 



1 + Ce 2 + [ C - y + 5] ei e 2 + (ho 2 - y + 1 



(63) 



At first order, this is nothing but the consistency condition of inflation, i.e. r = — 8n T . 
Therefore, any upper bound on the first slow-roll parameter e\ gives an upper bound 
on the contribution of gravitational waves. In fact, the ratio r itself is not observable 
because one must take into account the evolution of the transfer function when the 
modes re-enter the Hubble radius. If we formulate the ratio tensor to scalar in terms of 

CMB multipole moments, then one obtains 

r<(T) 



r 10 = -^ = -f 10 (h,n A ,---)n T . (64) 
°io 

For the concordance model, one has j\o ~ 5 which implies that T\q ~ 10ei. We see that 
an upper bound on e\ puts an upper bound on r 10 . This also implies an upper bound 
on the energy scale of inflation. Indeed, on very large scales [I <C 20), the multipole 
moments can be written as 

/7(B) ~ ZJ1 1 

and the measurement of the quadrupole C 2 by WMAP means that 

F V Q 0rcei ( Q *™- ps )\ (66) 



77T.pi J \ T 



where jHj 



Qrms-PS /5C2 in-6 

— ^— = ^^ 6 x 10 6 . (67) 

As announced, any upper bound on the first slow-roll parameter implies an upper bound 
on H/m Pl . In the next section, the WMAP third year data are used to constrain the 
slow-roll parameters and an upper limit on r 10 and H/m P1 is derived. 



3.2. WMAP data constraints on the slow-roll parameters 

In this section, before turning to the specific models of inflation we are interested in, we 
derive the bounds that the Hubble-flow parameters have to satisfy given the third year 
WMAP data (WMAP3). Following |S2], we used a modified version of the CAMB code [TO] 
to compute the CMB temperature and polarisation anisotropics seeded by the slow-roll 
scalar and tensor primordial power spectra of equation (|45|). The parameter space has 
been sampled by using Markov Chain Monte Carlo (MCMC) methods implemented 
in the C0SM0MC code JT] and using the likelihood estimator provided by the WMAP 
team IH El Uj- The likelihood code settings have been kept to their default values 
which include the pixel based analysis at large scales, a Gaussian likelihood for the beam, 
with diagonal covariance matrix, and point source corrections [HI ESI E] • We have also 
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checked the sensitivity of our results to modifications of these options and no significant 
deviations has been observed. The assumed cosmological model is a flat ACDM universe 
involving a minimal set of four cosmological base parameters: the number density of 
baryons Q^, of cold dark matter fldm, the reionization optical depth r and 9 which 
measures the ratio of the sound horizon to the angular diameter distance (see Ref. 
Moreover, we have only considered the WMAP3 data together with the Hubble Space 
Telescope (HST) constraint (H = 72±8km/s/Mpc [S3]) and a top hat prior on the age 
of the universe between 10 Gyrs and 20 Gyrs. The convergence of the chains has been 
assessed by using the Gelman and Rubin R-statistics implemented in C0SM0MC which is a 
measure of the variance of the means divided by the mean value of the variances between 
different chains EJ. Unless otherwise specified, the iterations have been stopped 
once R — 1 < 3%, which corresponds to a few hundred thousand samples depending on 
the model explored. 

As described in the previous section, the set of primordial parameters in the slow- 
roll approximation involves the value of the scalar power spectrum = k 3 P^ (k*) and 
the Hubble-flow parameters evaluated at the pivot scale k*. 

3.2.1. First order slow-roll expansion At first order in the slow-roll expansion, the 
relevant parameters are e± and €2 while the running of the spectral index vanishes. 
This accounts for three primordial parameters including the amplitude of the scalar 
power spectrum leading to seven model parameters in total. The best fit model has 
X 2 = 11252.2 (for comparison, the standard power law parametrisation almost leads to 
the same best fit value y 2 = 11252.4). 

In addition to the usual uniform top hat priors on the cosmological parameters ^1] , 
we have chosen a uniform prior on ln(10 10 Pj in the range [2.7,4.0], as well uniform 
prior on €2 in [—0.2, 0.2]. As discussed in [HE1EI!, a prior choice on e\ is not innocuous. 
This parameter encodes the amount of gravitational waves contributing to the CMB 
anisotropies and its order of magnitude is not known. Moreover, as shown in the next 
sections, even in the simplest single field inflationary models under scrutiny in this work, 
the orders of magnitude of ei may vary considerably. As a result, we have performed 
two MCMC analysis: one where e\ has an uniform prior in the range [10~ 5 , 10 _1 ] and 
the other with a Jeffreys' prior by choosing an uniform prior on log(ei) in [—5,0]. We 
have plotted in figure El the marginalised posterior probability distributions for the base 
cosmological parameters together with H and for convenience. The solid and dashed 
line correspond to the uniform and Jeffreys' prior on e%, respectively. 

The deviations induced by the prior choice on ei are the consequences of the 
degeneracy between the two slow- roll parameters (see figure EJ). As discussed in 
this is the result of both their influence on the spectral index n s = 1 — 2ei — €2 and 
the tensor mode contribution to the CMB anisotropies encoded in e\. Since a Jeffreys' 
prior on e\ gives more statistical weight to its small values, the accessible volume in the 
parameter space is enlarged in a region where e\ is small and the tensor modes have 
not observable effects. In some sense, this prior choice tends to favour pure scalar mode 
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Figure 2. Marginalised posterior probability distributions for the base ACDM 
cosmological parameters together with the cosmological constant and the Hubble 
parameter, obtained at first order in slow-roll expansion. The solid black lines 
correspond to an uniform prior choice on e\ while the dashed blue ones to an uniform 
prior on log(ei). The dotted black lines is the mean likelihood for the former prior. 



models. On the other hand, in the region where €\ has a significant effect, e 2 is pushed 
to lower values in such a way that the scalar spectral index remains compatible with the 
data while the cosmological parameters react to a significant tensor mode contribution. 

In figure we have plotted the marginalised posteriors for the primordial 
parameters. Similarly, the probability distribution on e 2 ends up being affected by the e\ 
prior choice. Since for small values of t\ the spectral index approaches n s ~ 1 — e 2 , under 
the Jeffreys' prior on e\ one finds the positive values of e 2 favoured: —0.02 < e 2 < 0.09 
at 2a level while —0.07 < e 2 < 0.07 under the uniform e\ prior. We also obtain the 
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0.03 



0.02 



0.01 




Figure 3. 68% and 95% confidence intervals of the two-dimensional marginalised 
posteriors in the slow-roll parameters plane, obtained at leading order in slow-roll 
expansion. The shading is the mean likelihood and the left plot is derived under an 
uniform prior on ei while the right panel corresponds to an uniform prior on log(ei). 




0.01 0.02 0.03 -5 -4 -3 -2 -1 

Figure 4. Marginalised posterior probability distributions for the primordial 
parameters in the first order slow-roll expansion. As in figure |2 the black solid curves 
are derived under an uniform prior on ei whereas the dashed blue curves correspond 
to a flat prior on log(ei). At two-sigma level of confidence one has the marginalised 
upper bound ei < 0.022. 
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marginalised upper bound e\ < 0.022 at 95% of confidence, or in terms of the observed 
tensor to scalar ratio and scale of inflation: 

no < 0.21 , < 1.3 x 10" 5 , (68) 

m P1 

again at 2a. 

Notice that without marginalising over 62 one obtains slightly weaker 2a limits 
e\ < 0.028 and r w < 0.28. These can be compared to the unmarginalised bounds 
coming from the WMAP first year data and given in reference [22], namely e\ < 0.032 
and rio < 0.32. 

3.2.2. Second order slow-roll expansion As discussed in sectionOJ at second order in the 
slow-roll approximation, one has to consider the third Hubble- flow parameter €■$. In the 
following, we reiterate the previous MCMC analysis on the parameters space enlarged 
by e 3 under an uniform prior choice in [—0.1, 0.1]. All the other priors have been kept as 
in the first order analysis, as well as the comparison between the uniform and Jeffreys' 
prior on e\. We are therefore dealing with a model involving eight parameters. 

As expected for well-constrained parameters, we find no significant deviation on the 
base cosmological parameters and scalar power spectrum amplitude between the first 
and second order slow-roll models. The two-dimensional marginalised posteriors for the 
Hubble-flow parameters are represented in figure El Similarly, the one and two-sigma 
confidence intervals in the plane (ei,e 2 ) are found to be slightly enlarged compared to 
their first order equivalents, but not more than what one may expect from the inclusion 
of a new parameter. 

The marginalised posteriors associated with the Hubble-flow parameters are plotted 
in Fig. El Apart slightly weaker constraints on e± and €2, one may notice a weak, but 
non-significant, tendency of running associated with positive values of £3. However, 
this effect appears only when a Jeffrey's prior is chosen on e\ and disappears under the 
uniform prior choice [521 El EH ■ 

3. 3. Large field models 

Large fields models (or chaotic inflation models although this is not very 
appropriate [SH]) are characterised by the monomial potential given by [UU] 

V ( ( j ) )=M 4 (^-Y . (69) 



ni 



pi 



There are only two free parameters, M and p and the energy scale M is uniquely 
determined by the amplitude of the CMB anisotropies, and thus the WMAP 
normalisation. This is probably the simplest inflationary scenario since, in the slow- 
roll approximation, everything can be integrated analytically. This family of potentials 
is represented in figure From the expression of the potential, the three slow-roll 
parameters can be easily computed and reads 
p 2 ml, p ml, 

ei= 16^T' £2= 4^' e3 = 62 - (?0) 
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Figure 5. One and two-sigma confidence levels of the two-dimensional marginalised 
posteriors in the Hubble-flow parameter planes (solid black contours). The dotted red 
contours correspond to the same confidence intervals derived at first order in slow-roll 
expansion (see figure OJ). The left panels are associated with an uniform prior on ei 
whereas the right ones correspond to an uniform prior on log(ei). As can be seen in the 
bottom right panel, positive values of 63 are slightly favoured, but in a prior dependent 
way while £3 = remains within the one-sigma contour. 



These slow roll parameters are represented in the two bottom panels in figure There 
are monotonic functions of <fi and they are decreasing as is increasing. One can 
immediately deduce that, for a given p, the model in the plane (ei, €2) is represented by 
the trajectory e\ = (p/4)e2. We can also estimate the value of the inflaton at the end 
of inflation, defined to be the time at which t\ = 1 (see before). This leads to 

0end P 



(71) 



Moreover, the slow-roll equation of motion leads to a solution which is completely 
explicit. Integrating the following quadrature 



N 



in 



ni 



V( X ) , 8tt 
dx = 



pi J $ 



<f>/m P1 



xdx . 



P J<j> in /m P1 



(72) 



Inflation after WMAP3 



20 





Figure 6. Marginalised posterior probability distributions for the Hubble- flow 
parameters up to second order in the slow-roll expansion. The solid thick black curves 
corresponds to a uniform prior on ei whereas the solid thin blue curves are derived 
under an uniform prior on log(ei). The dashed curves are the corresponding first order 
slow- roll posteriors of figure 0] 



leads to an explicit expression N = iV(0) which can be inverted and reads 

T>- (73) 




m P1 y \m Pl J An 
This expression also allows us to obtain the total number of e-folds. It is sufficient to 
evaluate the above expression for = en d, 0end being given by (J7TJ). One arrives at 

WT = W*.y p (74) 

p \m P1 J 4 

which can be very large if the initial energy density of the inflaton field is close to the 
Planck scale m pl . However, the model remains under control only if the initial energy 
density is smaller than m pl and this imposes a constraint on the initial value of the field, 
namely 

^~ImJ • (75) 

Let us notice that, when the inflaton energy density approaches the Planck energy 
density, quantum effects become important. In this case, the formalism of stochastic 
inflation must be used [HH E21 E31 EH ESI EEl EZl EH1 
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Figure 7. Top left panel: large field potential forp = 2 (solid green line), p — 4 (dotted 
blue line) and p = 6 (dashed red line). Top right panel: logarithm of the potentials 
for the same values of p (and the same line style code). The required flatness of the 
potential becomes obvious on this plot. Bottom left panel: slow-roll parameter e\ for 
a large field potential with p = 2. The shaded area indicates the breakdown of the 
slow- roll inflation (strictly speaking when the acceleration stops). Bottom right panel: 
slow- roll parameters £2 and £3 for a large field potential with p = 2. Only one curve 
appears because £2 = £3. On this plot, the shaded region signals the breakdown of the 
slow-roll approximation but not necessarily the end of the accelerated phase. 



We now turn to the explicit determination of the slow-roll parameters. We have 
seen that the model is represented by the trajectory e\ = (p/4)e 2 but observable models 
only lie in a limited portion of this straight line. Indeed, the Hubble-flow parameters 
should be evaluated when the scales of astrophysical interest today left the Hubble 
radius during inflation. Let us call the value of the inflaton at that time 0*. Then 0* 
can be expressed in terms of iV*, the number of e-folds between the time of Hubble exit 
and the end of inflation: 

/"(/Wd/mpi 

iV* = / xdx , (76) 

P J<t>*/mp\ 

from which one deduces 

4>t P fnr , P 



77lpi 47T V I 



iV* + 7 • (77) 




Figure 8. Slow-roll predictions for large field models in the plane (ei, 62). The dotted 
black lines indicate the lines of constant spectral indices from n s = 0.8 (on the right) 
to ii s = 1.1 (on the left), the red line corresponding to a scale invariant spectrum 
n s = 1. The solid blue lines separate the regions £2 < 0, £2 > 0, €2 < 2ei and 
£2 > 2ei associated with different energetic evolutions (see text). As expected, large 
field models lie in the region €2 < 2e\. The dotted-dashed blue curves indicate the la 
and 2(7 confidence intervals given the WMAP3 data (see section Each coloured 

segment represents the prediction of a model given p (see colour bar) and for a number 
of e-folds between the end of inflation and the Hubble exit varying in [40, 60]. 



Therefore, the slow roll parameters can be expressed as 

£l = 4(iV*+p/4)' e2 = i^T^' €3 = €2 - (78) 
The number of e-folds iV* between the end of inflation and the Hubble length exit can 
be thought as the arc length along the straight line representing the large field models. 
The value of N* can be calculated once one knows the entire history of the Universe, 
including the reheating phase [69] . In the following, we will consider that 

40 < JV, < 60 , (79) 

when dealing with the slow- roll models and in fact, as explained in 69J, one could 
consider an even more restricted range for the quartic model p — 4. 

The slow-roll predictions for the large field models are represented in figure [HI As 
expected, the whole family lies in the region > and 62 < 2ei. According to the 
previous discussion, this means that, as inflation proceeds, the kinetic energy increases 
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with respect to the total energy density and, at the same time, the absolute value of this 
kinetic energy density decreases. From figure |Hl all the models with p > 4 lie outside 
the 2a contour. The quadratic (or massive) model remains compatible with the data 
and predicts quite a high contribution of gravitational waves, up to r 10 ~ 10% level. 
Having found the compatible values of the parameter p, our next move is to estimate 
the numerical value of the parameter M. This can be done from the measurement of 
the CMB quadrupole fUJ) made by the WMAP satellite 



Qrms-PS _ 1 Hi 2 V* 

T 2 607re* m\ x 456* m% l 



JO) 



In the case of large fields model, this implies 



AfV_ (45/2)p y /z QL s -ps 



(AN, + p) p/2+1 V P J T 2 



and given the constraints on p and N*, this leads to 

4 x 10" 4 < — < 1.1 x 10~ 3 . (82) 
m P1 

We recover the conclusion that, for large field models, inflation take place close to the 
Grand Unified Theory (GUT) scale. 

To conclude this section, let us come back to the question raised after equation 
(I42J) . namely the error caused by expressing the slow- roll parameters e n in terms of the 
potential V and its derivatives. If we take into account the more accurate equation (|44|) . 
then the expression of the parameter ei for the quadratic model becomes 

ei = ^^{ 1 ~^^ + --)- (83) 
From this expression, if one tries to calculate a new value of e nd, one immediately faces 
the problem that the above expression leads to a second order algebraic equation with 
a negative discriminant, i.e. an equation which does not admit real solutions. This is 
not surprising since, near the end of inflation, the correction is a priori of order one. 
In order to have an order of magnitude estimate of the corrections, one can still work 
with the former value of en d in (|7ip. Then, as expected, one finds that the correction 
is of order 1/N 2 ~ 4 x 10~ 3 and, hence, does not in any way modify the conclusions 
established before about the compatibility of the large field models with the WMAP 
data. 



3.4- Small field models 



We now turn to another class of models, namely small field models. In this case, the 
inflaton potential can be written as [7U1 EH H2] 

pi 



V(<t>) = M 



1 - 



M) 



The potential is characterised by three parameters, the energy scale M, \x and the index 
p, i.e. one more parameter with respect to the large field models. 
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Figure 9. Top left panel: small field potentials for p = 2 (solid green line), p = 4 
(dotted blue line) and p = 6 (dashed red line). Top right panel: logarithm of the same 
potentials than in the top left panel. Bottom left panel: slow-roll parameter e\ for the 
model p — 2 with fj,/m P1 = 0.1 (solid line), /i/m P1 = 0.5 (dotted line) and fi/m P1 = 1 
(dashed line). The dashed area indicates where the slow-roll breaks down. Bottom 
right panel: slow-roll parameters ti (solid line) and 63 (dotted line) for fi/m rl = 1. 
Note that the slow-roll parameter 62 is non- vanishing in the limit (j>/fi — > where its 
value reads (m P1 //i) 2 p/(47r). 



The potential, as well as its logarithm, are represented in figure |H1 For these models, 
it is straightforward to calculate the first three slow-roll parameters. They read 



167T \ [i 



(MO 



2p-2 



£ 2 = — 



P_ 

4tt 



e 3 = 7T- 
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1 - (0//i)f ' 
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+ (p-l)(p + 4) 



+ (p-l)(p-2) 
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57) 



and they are represented in the two bottom panels of figure M As for large field models, 
they are monotonic functions of the field. However, these parameters are now increasing 
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as the vacuum expectation value (vev) of the inflaton is increasing during small field 
inflation. 

Then, we follow the same steps as in the previous subsection. Our first goal is 
therefore to determine when inflation stops. Requiring ei(0 en d) = 1 leads to an algebraic 
equation, namely 

endA 2p _ 2 ^y_^^y^y p - 2 +1 = _ (88) 



fl J \ fl J 167T \ fl J \ fl 

This equation can be solved explicitly for p = 2 and the solution reads 



Wd 1 m P1 ( 1+ I 16V (g9) 



jj, 4^/7? fi \ y m pl 



In fact, among the two solutions (j88jl admits, we have chosen the one with the minus 
sign. This is because, since the vev of the inflaton field is increasing as inflation proceeds, 
the accelerated phase stops for the "smallest" value of en d- Let us notice that, in the 
limit, fi/m Pl — > +00 one has en d - ► A 4 but we always have en d < fi, in other words the 
limit is approached by lower values^. On the other hand, if we now assume fi/m Pl <C 1, 
then the previous solution can be approximated by 



end ~ 2^-^- , (90) 



fi m P1 

and, therefore, under the assumption fi/m Pl <C 1 one has 4> e nd/ H «C 1. 

The above considerations were established for the case p = 2. If p 7^ 2 then, as 
already mentioned, equation (J85j) can no longer be explicitly solved but we can still 
approximate its solutions in the limits considered before. Indeed, for fi/m Pl ^> 1, 
equation (jSSjl reduces to y 2 — 2y + 1 = 0, where y = (0 e nd//-O p > the only solution of 
which is en d = /i. Consequently, as in the case p = 2, end — > jj, when /i/m P1 — > +00. 
On the other hand, since we always require en d/ V* < 1, the two terms (0 cn d/ n) 2p and 
(0end/Ai) P in (jHHj) can be neglected in comparison to 1, provided fi/m Pl 1. Then, 
keeping the third term on the left-hand side leads to 

1/O-2) 

(91) 



^end 



t> 



16n ( ji, 



p 2 \m ¥ 



and one can check that, for p = 2, this reproduces equation (|90|). If, on the contrary, 
one does not have fJ,/m Pl 1 (or /i/m P1 ^> 1), then equation (|SS|) can only be solved 
numerically. 

The next step is to obtain the classical field trajectory for these models. This can 
be done if the slow-roll approximation is satisfied but, even in this case, the classical 

I Of course, this limit is physically questionable because the expression ()84|l could viewed as a Taylor 
expansion in 4>/fi. As long as (j> ~ /a, the accurate form of the potential should be given and, 
consequently, other power of <^f [i considered |78| 
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trajectory can only be found implicitly. In terms of total number of e-folds N, we have 
for the small field potential 



giving 



N 
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for p 7^ 2, while for p = 2 one has 



N = 2vr 



mi 



ft 



(92) 



, (93) 



(94) 



For p — 2, one can invert the above relation and express in terms of the number of 
e-folds to get 



-w 



cxp 



2 iV / m, x 



27T \ /X 



(95) 



where Wo (a;) is the principal branch of the Lambert function j^j. This special function§ 
is the solution of the equation 

W{x)e w ^=x. (96) 

In our case, we have to choose the principal branch since </>///< 1. 

Finally, the last step consists in determining the link between iV* and *. For p ^ 2, 
using the classical trajectory obtained before, one has to solve the following equation 
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(97) 



where en d and iV* are known from the previous discussion. 

For p = 2, this equation can be explicitly solved and one gets 



\ 



-'end 



exp 



fend 



N* f m 



2n \ fi 



(9* 



In the limit fi/m Pl <C 1, the argument of the Lambert function is small and Wq(x) ~ x 
leading to 



-*cnd 



exp 



-*cnd 



2tt \ fi 



(99) 



On the other hand, in the limit /i/m P1 — > +oo then, as established before, en d — * p, 
and the argument of the Lambert function goes to — 1/e. Since W (— 1/e) = —1, one 
obtains 0* — > /i. 

§ It is also called ProductLog [ ] in Mathematica. 
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If p 7^ 2, then equation (J97|) cannot be solved explicitly. However, in the limit 
/i/m P1 <C 1, one has 0* /i and (fiend/ H "C 1 allowing (}9"T)) to be approximated as 

pN* ( m P1 \ 2 i 2 />cndV~ P ^ 100 ^ 



2-P V/ 1 
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2-p 



^end 



Using the expression of <fi eno \ obtained in (j91|) renders the second term in the right-hand 
side of the previous equation proportional to (m P1 / '^(p -2 )/^ -1 ). Under the considered 
limit yu/m P1 <C 1, this term can be neglected in comparison to the the first one which is 
proportional to (m P1 /yu) 2 . Putting everything together, one gets 

1/(2~P) 

(101) 
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Let us stress again that the above result is valid only for /x/m P1 <C 1. If this is not the 
case, one has to rely on numerical calculations to find the correct value of 

The values of the slow-roll parameters directly stem from the previous 
considerations. If p = 2, it is sufficient to use the value of 0* found in (|9*5j) into 
the expressions ()85j) to ([87)1 . The only underlying approximation being in that case the 
inversion mentioned in section l2~2l Besides, assuming the limit /i/m P1 <C 1, one arrives 
at 
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(103) 
(104) 



These expressions coincide with the formulas already used in the literature [721 113 IS] 
and it is crucial to keep in mind that they are valid for /i/m P1 <ti 1 only. 

In figure El the slow-roll parameters E\ and e<i stemming from the approximate 
equations ()102j) and (jK)Hj) are represented together with the ones obtained from the 
exact expression of 0* in (|98|). The plot is made for iV* = 40 (and p = 2). One sees 
that as long as /x/m P1 1, the two expressions of ej and e 2 are in good agreement but 
when [A ~ m PJ the difference is no longer negligible. It turns out that this difference is of 
utmost importance in view of the current data because using only the equations ()102|) 
to ()104j) leads to the incorrect conclusion that the small field model p = 2 is ruled out. 
In fact, using directly the slow- roll equations shows that this model is still compatible 
with the observations. Indeed, very roughly speaking, the WMAP3 data are compatible 
with e\ < 0.03 and e 2 — 0.05. If one decides to use ()103j) for then fi ~ 1.8m P1 and 
inserting back this value into ()102|) . one finds e x ~ 0.13, i.e. a value in tension with the 
observations. The true story is quite different. As can be seen in figure E3 the slow-roll 
value €2 — 0.05 is perfectly compatible with a value for ei satisfying the WMAP3 bound. 
The model is thus still compatible with the observational constraints. This conclusion 
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Figure 10. Slow-roll parameters e\ (blue lines) and £2 (red lines) for the small field 
model characterised by p = 2 and = 40, plotted for different values of the scale 
/j,/m P1 . The dotted lines are the approximations established in equations Ijl02|l to (|1()4|) 
whereas the solid lines are the slow-roll parameters expressed in terms of the Lambert 
function: the only approximation involved in that case being the inversion discussed 
in section l2~2l 



is indeed confirmed by the exact numerical integration performed in the next section. 
In addition, we will demonstrate below that /i/m Pl <C 1 implies a very low energy scale 
during inflation that can be, in some cases to be discussed in the following, below the 
nucleosynthesis scale, i.e. already ruled out. In this situation, the small field models 
make sense only if /i > m P1 , that is to say precisely the situation where it is necessary 
to carefully evaluate the values of the slow-roll parameters. 

Let us now turn to the expressions of the slow- roll parameters in the case p 7^ 2. If 
/i/m P1 1, then one has 
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(105) 



(106) 



N*p — 2 

which matches to the expressions usually used in the literature. On the other hand, if the 
limit /i/m P1 -C 1 is not satisfied, the slow-roll parameters can only be estimated from 
numerical calculations. In figure ^2 we have precisely compared the approximations 
given by the equations ()105j) and (jlOfij) to the directly computed slow-roll parameters 




Figure 11. Comparison of the slow-roll predictions for small field models with p = 3 
and 0.1 < n/rripi < 20 in the plane (ei, £2)- The solid red and blues lines, as well as the 
blue dotted-dashed contours, have the same meaning as in figure|Bl The approximated 
slow-roll values are obtained using the equations H105(l and (|1()6|) and are represented 
by the two red dotted vertical lines. The red dotted line on the left corresponds to 
models with different values of /i/m P1 at fixed TV* = 60 while the red dotted line on 
the right is obtained for =40. The correct slow- roll parameter values are obtained 
by numerical integrations and are represented by the line segments. Each line segment 
corresponds to a model with a fixed value of /z/m P1 , indicated by the colour bar, along 
which the quantity N* varies from 40 to 60. 



in a model with p = 3. The approximated values are represented by the two red dotted 
curves, the left one corresponding to the choice iV* = 60 and the right one to = 40. 
These two curves appear as vertical lines because the expression of e 2 in (jl05|) does not 
depend on /x/m P1 . The correct values for the slow-roll parameters obtained by numerical 
determination of 0* are represented by the line segments, iV* varying from 40 to 60 along 
each segment. Although both methods are in good agreement for yu/m P1 <C 1, they differ 
when /1 > m P1 . As it was the case for p = 2, we see that using only (jl05|) would lead 
us to the erroneous conclusion that the model p = 3 is compatible with the data only if 
/U / m P1 <C 1 . On the contrary, the correct values of the slow- roll parameters indicate that 
models with /i > m P1 are simply in perfect agreement with WMAP3 data. Therefore, 
given the accuracy of the current CMB data, it becomes mandatory to carefully estimate 
the slow-roll parameters in the case of small field models. This will be confirmed by the 
full numerical integration of the power spectrum performed in the next section. 
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Figure 12. Small field models in the plane (£1,62) for 0.2 < fi/m P1 < 20, p — 2 (top 
left panel), p — 4 (top right panel), p = 6 (bottom left) and p = 8 (bottom right). The 
solid red and blues lines, as well as the blue dotted-dashed contours have the same 
meaning as in figure [S] Each curved segment represents the slow-roll predictions along 
which varies from 40 to 60. 



Let us now more precisely consider the predictions associated with the small field 
models. In figure El the slow-roll parameters values are represented for 0.2 < /i/m P1 < 
20 and for four models, namely p = 2,4,6 and p = 8. As expected, the models are 
located in the region where €2 > 2ei > 0, i.e. in the region where the kinetic energy grows 
with time during inflation. One can see in figure that, for any value of p (including 
p = 2), the models are in good agreement with the data if fi > m P1 . This conclusion is 
confirmed by the full numerical integration done in the next section. In fact, only the 
very small values of ^/m P1 are problematic since they end up being associated with a 
too large slow- roll parameter €2- This can be seen in the top left panel (p = 2) where 
the line segments representing the /i/m P1 < 1 models are only marginally compatible 
with the confidence intervals or even ruled out. In fact, in the top left panel, this is not 
so apparent at first sight since the models corresponding to the segment line which lies 
between the one and two sigma contours are still acceptable fits to the data. However, 
one has to realize that, for even smaller values of yu/m P1 , the corresponding line segments 
are actually outside the figure, hence the above claim. Moreover, one notices that the 
segment line has the tendency to become vertical which indicates that the value of €2 is 
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independent from iV*. This is in full agreement with equation (J103|) giving e<i- 

Finally, one word is in order on the CMB normalisation in those models and the 
values of e±. It is known that the very small values of €\ imply a very small contribution 
of gravitational waves. But they also imply a quite low energy scale during inflation. 
This allows us to derive a lower bound on e%. Indeed, the normalisation of the spectrum 
is given by 

V * ~, 4561 ^rms-PS nn? N 

m pl " 2 T 2 ' 1 ' 

Now, it is physically motivated that inflation must take place at an energy scale at 
least greater than, say, the TeV scale. This means V^/m^ > 10~ 64 and for p = 2 (and 
AT* = 50), one gets /i > 0.25m P1 . Such a strong lower bound on /i is essentially due to the 
exponential behaviour of the slow-roll parameter with respect to the model parameters. 
In fact, this condition can equally be worked out for p ^ 2. For instance, with p — 2.1 
one obtains fi > 0.03m P1 , p = 2.5 leads to /i > 7 x lCT 6 m P1 , p = 3to/i>3x 10~ 9 m P1 and 
p = 10 corresponds to fi > 8 x 10 _21 m P1 . These results sustain the remark made before 
motivating a careful determination of the slow-roll parameters in the regime /i > m P1 . 



3.5. Hybrid inflation 

Let us now turn to hybrid inflation. This case is slightly different since hybrid inflation 
is in fact a two-field model with the potential I7H| ITT)] 

V (0, V) = \m 2 <P 2 + j - A 2 ) 2 + , (108) 

where is the inflaton, ip the waterfall field and A' and A are two coupling constants. 
Inflation proceeds along the valley given by if) = and, in this case, the potential reduces 
to an effective single field potential that can be written as 

pi 



V{<t>) = M 



1 + 



(109) 



with p = 2 and where one has used the following redefinitions 

A /1/4 A /A 7 A 2 . . 

M ^^T' " = V2m- (110) 

As for small field models, the effective potential (jl09|) depends on three parameters, 
namely M, fi and p. In fact, as mentioned before, p = 2 for the two field model given in 
()108|) but, to consider the most general situation we leave p unspecified in the following 
equations. Indeed, one could consider a model where the inflationary valley is described 
by, say, a quartic potential since the instability mechanism is independent from the 
shape of the valley. Moreover, in the second part of this article, when we perform the full 
numerical analysis, we will consider p as a free parameter and obtain its corresponding 
posterior probability distribution. Finally, let us emphasize that since we use a single 
field modelisation of hybrid inflation, such a representation cannot account for multifield 
effects as the generation of isocurvature modes or cosmic strings [77J [7HJ [T71 [7^]. The 
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Figure 13. Top left panel: hybrid field potential for p — 2 (solid green line), p = 4 
(blue dotted line) and p — 6 (red dashed line). Top right panel: logarithm of the 
potential for p — 2, 4, 6 (same colour code as in the top left panel). Bottom left panel: 
slow-roll parameter t\ for p — 2 and /i/m P1 = 0.1 (solid line), fi/m Pl — 1/(4a/7t) ~ 0.14 
(dotted line) and /x/m P1 = 0.3 (dashed line). Bottom right panel: slow-roll parameter 
£2 (solid line) and €3 (dotted line) both for fi/m P1 = 0.3. 



potential f)109j) and its logarithm are represented in figure for different values of the 
power p. 

The calculation of the slow-roll parameters proceeds as usual. The explicit 
expressions are given by 
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The three parameters are plotted in figure in the case p = 2. We see that the 
parameter e\ has a maximum at <p/fi = 1 which corresponds to the inflexion point of 
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In V . If fi < 1/ (4-^/7r) then e\ < 1 for all values of 0. We come back to this point below 
when we discuss how to stop inflation. Another specific feature of hybrid inflation in 
comparison to large and small field models is that, as can be seen on the bottom right 
panel, the parameters e 2 and e 3 can be negative. In particular 

hm £2 = _ KP-1) (^\ 2 (iY-\ (1U) 



0/M->O 47T \ yU / 

and e 3 blows up in the limit (4>/ f) p ~ > P ~ 1- 

We now discuss how inflation ends and how to calculate e nd- In the hybrid 

scenario there are a priori two mechanisms for ending inflation. Either inflation stops 
by instability when the inflaton reaches a value 

cri = yA, (115) 

for which the mass in the direction perpendicular to the inflationary valley becomes 
negative, or the slow- roll conditions are violated and e± = 1. The latter condition 
happens for a field value (fi e solution of 

The difference between this equation and (JHHj) resides only in the sign of the (4> e / fj,) p 
term. This equation admits a solution only if the following condition is fulfilled: 
/i/m P1 < l/(4 v /7r). In figure IT3*l such situations correspond to the case where ei can be 
greater than one. If fi/m Pl > l/(4^/7r) then e\ is always smaller than one and inflation 
can only stop by instability. The above equation (jll6|) cannot be solved explicitly unless 
p = 2. In this case, one obtains 



mpl 'l±Jl-^|. (117) 



Of course, one recovers the fact that the solutions exist only if fi/m Pl < 1/(4-^). The 
positive sign corresponds to the largest root, while the minus sign corresponds to 
the smallest one, <fi~ (see figure fT3|) . In the limit \xjm vx <^ 1 the previous equation takes 



the form 

4>e m Pl 



'1181 



and one recovers the large field value of en( j [see for instance equation ([7Tjl]. More 
generally, in the limit of small /i/m P1 , one can even approximately solve (|116|) for p ^ 2. 
However, contrary to the small field models case and besides the term proportional to 
p 2 , one should keep the term (0 e //i) 2p rather than 1. This leads to 

— —- A 7=> ( 119 ) 

which is the large field model expression of </> en d for a potential with an arbitrary power 
p in the field. Therefore, the final value enc j of the inflaton is the maximum of cr i and 
(j) e . To decide which mechanism is realised in practice requires the knowledge of the 
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model parameters. However, one crucial interest of hybrid inflation is that inflation can 
proceed for small values of the inflaton vev. As we have seen before, if the inflaton vev is 
large, then the model is equivalent to a large field model which was already considered in 
a previous subsection. Therefore, in the following, we will be focused on hybrid inflation 
taking place for <fi < 4>~ only: the so-called vacuum dominated regime. Since €\ is always 
lower than unity in such cases (see figure H~3*|) . hybrid inflation must stop by instability. 
As a result, cr i < <fi~ will be considered in the following as a free parameter accounting 
for a total of four inflationary parameters. 

The next step consists in calculating the classical trajectory. Straightforward 
manipulations similar to the ones performed for small field models, lead to 
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for p y£ 2, where k = v87r/m P1 is the reduced Planck mass. For p = 2 one has 
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which can be inverted to express <fr i n terms of the number of e-folds: 
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(122) 



Once again, Wq(x) denotes the principal branch of the Lambert function [73]. The 
above expression is very similar to equation . except for the signs. 

As for the other models, the last step consists in determining the link between 
and 0*. For p ^ 2, using the classical trajectory obtained before, one has to solve the 
following equation 
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where cr i and N* are known from the previous steps. In general, this equation can only 
be solved numerically. However, if p = 2, one gets 
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The slow-roll parameters for p = 2 are then obtained by inserting the above equation 
into the expressions (jlll|) to ()113|) . We therefore have an explicit form for the three 
parameters €\, e<i and e% which become explicit functions of /x, cr i and iV*. In the so- 
called vacuum dominated regime we are interested in, the term (0//x) 2 in the expression 
of the potential (|1U9|) ends up being a small correction only. Note that in the opposite 
situation, i.e. the inflaton dominated regime, the model is equivalent to chaotic inflation 
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which was already treated in the section devoted to large fields models. This question 
is discussed in more details below. 

In Fig. we have plotted the slow-roll predictions in the plane (ei, €2) for different 
values of /i/m P1 and 10~ 6 < (fi CT i/ fi < 1CT 1 . The values of fi/m Pl that we consider go 
from fi/m P1 = 0.6 (top left panel) to fi/m Pi = 1.8 (bottom right panel). Let us now 
discuss in more details these plots. 

For fi/m Pl = 0.6, some models lie in the region €2 > and others in the region 
£2 < 0. The first ones corresponds to (p CT \/ ' \i close to the upper limit chosen for the plot, 
namely 10 _1 while the second ones are valid for small values as one can check on the 
colour bar. The interpretation is as follows. When the value of cf) cr \ is close to its upper 
bound, inflation proceeds in a regime where the term ((fi/n) 2 in the potential (jl09|) is 
dominant. In this case, the model behaves as a quadratic large field model and this 
explains why the models are concentrated along the line £2 = 2e L . When the value 
of (fieri/ fi becomes smaller, inflation can proceeds for very small values of the inflaton 
vev and one enters the vacuum dominated regime which is the main characteristic of 
hybrid inflation. Consequently, as can be seen in figure ITU the points representing the 
corresponding models leave the large field region €2 < 2t\ and penetrate in the region 
£2 < 0. In this region, the kinetic energy decreases absolutely but also relatively with 
the total energy density. With /i/m P1 = 0.6, these models are now disfavored by the 
WMAP3 data at more than 2a. 

However, the situation can be different when one starts to increase the value of 
n/m Pl . As we see in the top right panel of figure El arid even more clearly in the 
remaining panels, the models are no longer equivalent to quadratic large field models 
and the line €2 = 2e\ is never reached. This effect is amplified as fi/m Pl increases 
since then, the relative importance of the term ((fi/fj,) 2 is diminished. These remarks are 
nevertheless valid provided the range of (fi CTi is kept fixed. Clearly, as fi/m Pl increases, one 
could always increase the upper limit on cr i/ \x to counter-balance the above-mentioned 
effect. As can be seen in the left and right middle panels, all the models are almost 
concentrated in the vacuum dominated region €2 < and have the tendency to produce 
a blue tilted power spectrum, a standard characteristic of hybrid inflation. Some models 
with t2 < are still perfectly compatible with the data as can be seen in the middle right 
panel. Then, as /x/m P1 is further increased, the models become excluded due to their 
too high blue spectral index (bottom left and bottom right panel). Again, in this case, 
to re-obtain models in agreement with the data, it would be necessary to modify the 
upper bound on (fi CT i/fi. The parameter space of hybrid inflation and its compatibility 
with the CMB data is directly explored in the following section where the power spectra 
are numerically integrated. 
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Figure 14. Top left panel: slow-roll predictions for hybrid inflation with p = 2, 
/z/m P1 = 0.6, 10~ 6 < </>cri/M < 10 -1 . The colour bar indicates the value of log(0 cr i//i) 
and the number of e- folds is between 40 and 60. The other panels display different 
values of /i/m P1 , the other parameters being unchanged: fi/m P1 — 0.8 (top right), 
/Lt/m P1 = 1 (middle left), (J,/m P1 = 1.4 (middle right), [ijm vl = 1.6 (bottom left) 
and /x/m P1 = 1.8 (bottom right). The solid red and blues lines, as well as the blue 
dotted-dashed contours, have the same meaning as in figure |H1 
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3.6. Running-mass inflation 



The last type of model that we consider is the running-mass model (RM) |80 | l8T j 182]. 
Using the same parametrisation as before, the potential reads 



V(<j>) = M 4 
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■- + ln- 
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which is a function of 0//i only. This potential is characterised by three free parameters, 
M, /i and A. However it turns out to be more convenient to use a slightly different 
parametrisation (essentially to facilitate the comparison with the existing literature). 
For this reason, we re-write A as A//i 2 = k 2 c/2. Then, the potential takes the form 



V{<j>) = M 4 
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2 yj 0/ 

where 0o = /•*• 111 this expression, M, c and 0o are free parameters. Let us recall that 
c can be positive or negative |H] while = 0o is an extremum of V(4>), a maximum 
if c > and a minimum if c < 0. The potential and its logarithm are represented in 
figure fTol 

Running mass inflation can be realised in four different ways [80J, denoted as RM1 
to RM4, according to where the vev of the inflaton field is located along the potential, 
see figure ITBI (top panels). RM1 corresponds to the case where c > and < O . In 
this case, decreases during inflation. RM2 also corresponds to c > but with > O 
and increases during inflation. RM3 refers to the situation where c < and < 0o all 
the time. In this case, increases during inflation. Finally, RM4 has c < and > O 
decreases as inflation proceeds. Using the potential ()126|) . one can calculate the three 
slow-roll parameters e\, e<i and t^. Their explicit expression read 
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The slow- roll parameters are represented in the bottom panels in figure 

Let us now study how inflation stops in these models. A priori, the end of inflation 
is found from the condition ei = 1 and for this reason, it is interesting to look at 
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Figure 15. Top left panel: running mass potentials for Ktfio — 1 and c = 0.2 (dotted 
green line) and c = —0.2 (dotted blue line). The vertical dashed line indicates the 
position of </>o . Top right panel: logarithm of the potential (same colour code) . Bottom 
left panel: slow-roll parameter e\ for n<po = 0.1 and c = 0.6 (solid green line), c = 0.2 
(dotted green line), c — —0.6 (solid blue line) and c = —0.2 (dotted blue line) with 
K(f>Q =0.1. Bottom right panel: slow-roll parameters £2 and £3 for c = 0.6 (blue lines, 
solid for £2, dotted for 63) and for c = —0.6 (green lines, solid for £2, dotted for £3) 
with k0o = 0.1. 



the behaviour of e% in more details (see the bottom left panel in figure IT5|) . First of 
all, one notices that e\ — for (f> = O (at 0.1/ k in this plot) which marks the limit 
between RM1 and RM3 on one side and RM2 and RM4 on the other hand. Secondly, 
we also remark that the curves corresponding to a model with the same value of |c| 
are almost identical when ^ < 1, This comes from equation (| 12 7)1 where, in this 
limit, the denominator approaches unity while the numerator depends on c 2 only. If we 
approximate the denominator by one, then the maximum of €\ in the range < <fi < </>o 
is located at <fi ~ (fio/e for which e\ ~ c 2 k 2 0q/(2c 2 ) <^ 1, assuming the physical values 
of the parameters c and 0o we are interested in. This means that, for RM1 and 
RM3, one always has ei < 1 and inflation cannot stop by violation of the slow-roll 
conditions. Therefore, one must use another mechanism and, naturally, we will consider 
that inflation ends by instability at some critical value cr i. 
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In the regime where > 0o, corresponding to RM2 and RM4, the curves 
representing e\ for different c but the same |c| separate. From ()127|) . one sees that 
this is due to the influence of the denominator. In the bottom left panel in figure ITH| 
one notices that, for RM4, the situation is very similar to what was discussed before, 
i.e. inflation cannot stop due to lower than unity €\ values. Hence, one must also use 
the instability mechanism for this model. 

Finally, it remains RM2. For this model, the inequality ei < 1 can be violated 
and inflation could stop normally. In practice, this happens for large values of k0 and 
meanwhile the other slow-roll parameters have already reached values greater than one 
meaning that the slow-roll approximation has already broken down. However, as long as 
ei < 1, inflation is still proceeding. In the following, we will also assume that inflation 
stops by instability (inverted hybrid mechanism). 

In summary, for the four running mass models, we always consider that inflation 
ends at some value cr ; which is therefore viewed as an additional free parameter. Let 
us notice that, even if ei <C 1 while inflation is proceeding, one could have 62 > 1 at 
some point [80J. In such a case, inflation would not stop but the slow-roll approximation 
would break down when €2 = 1. As a result, the running mass inflation models under 
scrutiny account for a total of four primordial parameters. 

We now turn to the calculation of the slow-roll parameters. Our first step is to 
obtain the classical trajectory, that is to say the number of e-folds in terms of the vev 
of the inflaton field. In the case of the running-mass model ()126|) . it reads 
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where the exponential integral function [2U E2] is defined by Ei(x) = — &t€~ t /t. 
This expression cannot be inverted explicitly. However, in the limit k0 <C 1 (which is 
necessary to theoretically justify the shape of the running mass potential), the above 
expression can be approximated by 

(131) 

This form allows an explicit expression of the inflaton vev as a function of the number 
of e-folds, namely 
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From this expression, it is straightforward to calculate 0*. Remembering that inflation 
is supposed to stop by instability, one arrives at 
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Equipped with this value, it is then sufficient to evaluate the slow-roll parameters for 
this value of the inflaton vev using the equations (|127|) . (|128|) and (|129|) . 
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Figure 16. Top left panel: RM1 model with c = 0.01, k</> = 0.001 and 10~ 9 < 
rc^cri < 10~ 31 . Top right panel: RM2 model with c = 0.01, k</> = 0.001 and 

1-2.9 



io- 



< /t^cri < 10 -1 . Bottom left panel: RM3 model with c = -0.01, K<j) = 0.001 
and 10~ 5 



< K(f> cri < 10 -3 ' 1 . Bottom right panel: RM4 model with c = -0.01, 
Kifio — 0.001 and 10 -2 9 < rt<^ C ri < 10 _1 . The solid red and blues lines, as well as 
the blue dotted-dashed contours, have the same meaning as in figure [5] 



It may be interesting to present some analytical estimates of the slow-roll 
parameters. Instead of using the rather long equations (jl27j) . ()128j) and ()129|) . one 
can approximate them by assuming that the denominator is just given by M 4 , as done 
in 182j. In that case, the two first order slow- roll parameters read 
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In fact, it turns out to be convenient to define the following quantity 



s = cln 



(134) 
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which can also be written as 

s = -ce cN *ln(^j . (137) 

For RM1 and RM4, s > while for RM2 and RM3 one has s < 0. This quantity can also 
be used to estimate 4> cv \ since 4> CT i ~ 0oexp[— sexp(— a/V*)/c]. Then, the approximate 
equations giving the first two slow-roll parameters can be re-written as 

6i ^ y (^o) 2 e- 2fi / c , e 2 ~ 2c (l - £ + ^e~ 2 ^ . (138) 

The last equations means that the trajectory in the plane (61,62) can be expressed as 
£2 — 2(c — s) + 4ci/(k0o) 2 - If w e neglect e\ (see below), then one recovers the formula 
already derived in jHHl IHT| I82j. namely n s — 1 ~ 2(s — c). The same route for the third 
slow-roll parameter gives 
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(139) 



and if, again, one neglects ei, making use of equations and (|52*|) gives the the scalar 
running a s ~ 2sc. 

In figure we have represented the slow-roll predictions for the four versions of 
the running mass model in the plane (€1,62), together with the la and 2a WMAP3 
confidence intervals. The top left panel corresponds to the RM1 model with c = 0.01, 
K(f) = 0.001 and 10~ 9 < K<p cri < 10 -31 . The e\ parameter appears to be extremely small 
and, hence, the spectral index is approximately n s — 1 ~ 6%. According to the value of 
ft</> cr i, the spectral index can either be red for "large" values of 4> CV i relatively 

close to K(p , or blue for "small" values of K0 cr i. The relation n s — 1 ~ 2(s — c) also reads 



n s — 1 ~ —2c 



e cN * In I ^ ) + 1 



(140) 



and allows us to understand the behaviour of the spectral index. 

For the RM1 model, <p CT \ < <fio and the logarithmic term in the above equation 
is negative. When cri < O > this term is small, the constant term dominates in the 
squared bracket and, since c > 0, the spectral index is red. When cri O , the 
logarithm dominates, the squared bracket is negative and the spectral index becomes 
blue. 

The top right panel corresponds to the RM2 model with c = 0.01, k0o = 0.001 
and 10~ 2 ' 9 < K0 cri < 10 _1 : the spectral index is always red. This can be interpreted by 
means of equation ()140|) . Indeed, we now have cr j > (fio and hence the squared bracket 
is always positive ensuring that the spectral index remains lower than unity. Moreover, 
the larger cri , the redder n s in agreement with what is observed. 

The bottom left panel represents the RM3 model with c = —0.01, k0o — 0.001 and 
10~ 5 < K(f) C ri < 10 -31 . The spectral index can be red or blue depending on the value 
of K(f> cri . Since cr j < 0o the logarithm term is negative. If it dominates (0 cr j 0o) 
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then the bracket in equation (J14U)) is negative and therefore the spectrum is red since 
— c > 0. On the contrary, if the constant term dominates (0 cr j < 0o), then the spectrum 
can be blue. 

Finally, the bottom right panel corresponds to the RM4 model with c = —0.01, 
K(f>Q = 0.001 and 10~ 2 ' 9 < K0 cri < 10 _1 . Since (p CT \ > 4>o and — c > 0, the spectral index 
is always greater than one. This model can be compatible with the data only for "small" 
values of (fi cri as can be seen in figure HH1 

The previous observations agree with the existing literature, in particular with 
reference |80j . 

4. Testing exactly the inflationary models 

In this section, we do not use the slow-roll approximation but integrate numerically both 
the background evolution and the cosmological perturbations. This approach allows an 
exact determination of the power spectra of scalar and tensor modes assuming only 
the linear perturbation theory in General Relativity (BUI IHH EH1 EE1 E|- For a given 
model, these power spectra depend on the parameters characterising the potential and 
introduced in the previous sections. By coupling this mode by mode integration during 
inflation to a modified version of CAMB , we can use the MCMC techniques implemented 
in C0SM0MC to derive the constraints these parameters have to satisfy given the third 
year WMAP data. Let us stress that this method allows us to get marginalised posterior 
distributions directly on the potential parameters, out of any intermediate assumption. 
This approach is therefore different to the ones used so far and includes by construction 
a marginalisation on the reheating. As a result, the constraints we obtain include the 
effects coming from varying the number of efolds at which the cosmological perturbations 
can be generated [HHl E2 EH El El CTJ ■ 

The base cosmological parameters and their priors, as well as the HST data and 
the top age prior, are already described in section PT~2l Therefore, in the next section, 
we only discuss the method used to sample the inflationary parameters. 

4.1. Method 

Let us sketch how the numerical integration is performed for a potential of the form 

V{4>) = M 4 U((f)). 

The first step is to integrate the background evolution and this is done using the 
number of e-folds N as time variable. The energy scale M is initially set to an arbitrary 
non-physical value M = 1 (see below). The initial conditions 0i n and d0j n /diV are 
chosen such that there are at least 60 e-folds of inflation, the end of inflation being 
defined to be the time at which the exact Hubble flow parameter t\ — 1. Provided this 
condition is fulfilled, the initial conditions are in fact irrelevant thanks to the presence of 
the inflationary attractor [IT] . Once the background integration performed, the function 
4>(N) is numerically known in the range N G [0, N T ] (where N T > 60 is the total number 
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of e-folds during inflation) as well as all the other background functions, like the Hubble 
parameter H(N). 

In a second step, the equation (jZHj) controlling the evolution of the perturbations 
is numerically solved. This equation is fully determined only once the time-dependent 
frequencies u SjT (k, N), given in (|3()j). are known. This requires the knowledge of H(N) 
and its derivatives (up to third order), which simply comes from the background 
integration discussed above. In addition one needs the choice of some comoving 
wavenumbers k that will be the ones of astrophysical interest today. As a result, it 
is compulsory to be able to relate a comoving scale k during inflation to a physical scale 
fcphys defined at the present time. For that purpose, the complete history of the Universe 
needs to be specified and we now describe in more details how this can be achieved. Let 
us notice that equation (J29|) . written with N as time variable, takes the form 



dVs.T 1 d7Yd/i S)T 
diV 2 HdN dN 



k_^ 2 

n 



/i S ,T = 0, (141) 



where V M (iV) is the effective potential for the cosmological perturbations. We see that, 
besides the background evolution, we need to know the quantity k/TC during inflation 
as a function of the time variable N. Let us assume that we are given a physical scale 
today, say fc*/ao (in Mpc~ x for instance), where do is the present day scale factor. Then, 
one has 

k* kit do dend T N —N /-i a™ \ 
= ; r X X — = — e , (142) 

H a H(N) a cnd a(N) H(N) K J 

where a en d is the scale factor at the end of inflation and 

T = ^^, (143) 

a a end 

a constant which depends only on A;* and on the history of the Universe through the 
ratio ao/a en d- In order to evaluate this constant we assume that, after inflation, there is 
a period of reheating as described in section 12.1.21 characterised by the two parameters 
w re h and iV re h, followed by a radiation dominated era that can be described by fi° ad , 
followed by the matter dominated era. Let us be more precise about the reheating 
phase. It is interesting to consider the quantity i? ra d defined by 

lnR rad = In I J - 7 In ( J = - (-1 + 3w rch ) N ieh 



arch J 4 VPend/ 4 

_1 / 3 + 3 Wreh \ 

4 \5-3w reh J ' 1 ; 

where a re h and p re h are respectively the scale factor and the energy density in radiation 
at the end of the reheating phase. This last quantity is obtained from equation (f2U|) 
evaluated at t rch = r -1 , the last term in the squared bracket being neglected. From 
the previous definition, one notices that lni? ra d is exactly zero if w TerL = 1/3. In this 
case, the reheating phase cannot be distinguished from the radiation dominated era. If 
— 1/3 < U7 re h < 1/3 then lni? ra d < while if 1/3 < w rcrL < 1 then lni? ra d > 0. For 
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reasons that will be explained below in more details, it turns out to be convenient to 
define another quantity R 

In R = In R rad + ~ In (Vpend) . ( 145) 

The quantity p en d, the energy density stored in the scalar field at the end of inflation 
(ei = 1), is also completely determined once the background evolution has been solved. 
Therefore, if lni? ra( j is known, one can deduce \nR and vice- versa. Then, the ratio 
ao/dend can be written as [EH] 



a en d ti 
where p° ri is the critical energy density today, i.e. 

kVch = 3k 2 H 2 . (147) 

Therefore, given p° ri , f2° ad and R, i.e. a model of Universe between the end of inflation 
and the present time, one can calculate the constant T and hence all the terms in 
equation (j!41)l are explicitly known. In practice, we implement InR as the new 
inflationary parameter associated with the reheating and that is sampled from the 
MCMC: an uniform prior has been assumed on hiR in the range [lni? m i n , lni? max ]. 

Let us discuss how the limits -R m i n and -R max are chosen. First of all there are limits 
on p end . In order not to spoil the success of the Big Bang Nucleosynthesis (BBN) it is 
reasonable to require p en d > p nuc - Roughly speaking, this means p en d > 10 _85 m pl . On 
the other hand, one must have p en d < m ti i n order for the whole theoretical framework 
to be meaningful. In practice however, we do not need to implement this upper bound 
because the constraint H[ n f/m Pl < 1.3 x 10~ 5 derived in the slow-roll section (and coming 
from ei < 0.022) shows that the only viable perturbations have to verify p end < 10~ 10 m pl . 
In other words, we have 

- 187 < In (Vpend) < -20 , (148) 

the upper bound being not required a priori but being a consistency check that we 
should recover for the viable inflationary models. Now, if we use our reheating model, 
-R ra d can be rewritten as 

lni? rad = \~ 3Wreh > ( ^) - — -i -In ( 3 -±^I±) , (149) 

12(l + w reh ) VPend/ 3(l + w rch ) \5-3w rch J 

where p nuc < p reri < p en d and —1/3 < u> reri < 1 in order to satisfy the strong and 
dominant energy conditions. In this range, it is easy to see that the minimum is obtain 
for w re h = —1/3 and p re h = Pnuc whereas the maximum for w reri = 1 and p reri = p nuc - 
Finally, if we use the link between i? ra d and R, one arrives at 
1 11 

- In (Vpnuc) < lni? < -— In (fi; 4 pnuc) + 3 ln ( K Vnd) , (150) 

up to negligible factors depending on u; ren only. In the above equation, the bounds are 
our definition of lni? min and lni? max . Notice that these bounds explicitly depends on 
Pend which varies with the inflaton potential parameters. As a result, a "hard prior" 
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has been coded in COSMOMC to reject any sample involving a InP value that would not 
satisfy equation (jl50J) . Note that from the knowledge of InP and p en d, one can also 
derive lnP ra d and ao/a en d- 

From a chosen set of these parameters, the numerical integration can be 
straightforwardly performed to get the amplitude of the scalar and tensor power spectra 
at the end of inflation for any comoving wavenumber k that can now be related to their 
corresponding physical scales today 

However, as mentioned at the beginning, this result has been obtained with an 
arbitrarily chosen potential scale M = 1. In fact, the MCMC exploration is performed 
by directly sampling on the scalar power spectra amplitude P* at a given scale k* rather 
than on M directly An uniform prior has been chosen on ln(10 10 P !(! ) in the usual range 
[2.7, 4]. However, one has to restore the consistency between the value of P* and M since 
the former is uniquely determined by the latter. As we show in the following, a simple 
rescaling of the relevant functions can be used for this purpose. This rescaling has the 
advantage of being analytical and exact in the framework of the linear perturbation 
theory. 

Let us consider the following rescaling for the scale M or, equivalently, for the 
potential 

V{cj>) - sV{<f>) , (151) 

where s is a constant. What are the consequences of this rescaling on the other 
quantities? From equation ()11|) . one sees that the field <p(N) and its derivative are 
unaffected because only the logarithm of the potential appears in this formula. On the 
other hand, if we write the Friedman equation as 

"V, (152) 



H 2 



1 



K 2 



dN 



3 

then one notices that the Hubble parameter transforms as H — > s l ^ 2 H which 
immediately implies that p e nd — > sp en d- This also means that 7i = aH changes as 

n -> s^n. 

The previous considerations are valid for the background. Let us now see how the 
perturbed quantities are affected. The effective potential for cosmological perturbations, 
Vfj, in equation (|141|) . is invariant under the rescaling (|151|) . This can be seen for instance 
on the gravitational waves where V^(N) = a" /a = 1 — 7i ! /7i 2 . Therefore, requiring 

k-^s 1/2 k, (153) 

is sufficient to render equation (j!41j) invariant. However, this does not mean that the 
amplitude itself /i(fc, rj) is not changed. Indeed, although the equation of motion is 
invariant, the initial conditions are modified since they read 



^nrH(-) . (154) 



Because we deal with a linear equation, this implies that 

H {k, N) -> s-^fj, (s 1/2 k, N) . (155) 
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Therefore, the scalar power spectrum evaluated at the end of inflation, i.e. at N = N T , 
transforms in a very simple way, namely 

fi{k,N T ) 2 , gk3 ti(s 1 /%N T ) 12 



k 



(156) 



The last step is to determine the number s required to restore the consistency between 
our initial arbitrary normalisation M = 1 and the values of P* singled out during the 
MCMC exploration. From the previous discussing, the required value of s is simply 
given by the ratio P*/P^ M=1) where p( M =v is the amplitude of the scalar power spectrum 
stemming from the mode by mode numerical integration with M — 1 and evaluated at 
k = k^s" 1 / 2 . Then, the final power spectra do not need to be recomputed from scratch 
but can be directly deduced from the transformation (|156j) . and a similar one for the 
tensor modes. 

In principle, the previous procedure is well-defined but the fact that the wavenumber 
k must be continuously rescaled makes it quite difficult to implement in practice. In 
fact, it turns out to be more convenient to use a last trick that we now describe. The 
quantity p rc h can be viewed as an independent model parameter and, therefore, we also 
have the freedom to consider another model corresponding to a rescaled p re h- Let us 
consider the following transformation 

Prch -» S 2 preh , (157) 

where s is precisely the scaling number used in the above discussion on the rescaling of 
M. From these relations, one can immediately check that In R (or simply R) is invariant 
whereas P ra( j — > s -1 ' 4 Prad- From equation ()146|) . this also implies that 

a*_^ a ipa*_ t (15g) 

^end ^end 

If we now consider the quantity T we see that its transformation under the previous 
rescaling is the same if we change the scale k and leave p re h fixed or fix the scale k and 
change p re h- This is the main reason why it is more convenient to perform a MCMC 
sampling on the parameter R rather than P ra d since at a given R, the rescaling required 
on M does not induce a rescaling of the wave numbers. 

In the following, we have used such an exact numerical integration of the 
cosmological perturbations to constrain the models previously discussed in section 01 
This leads to new insights that we now describe and interpret. 



4-2. Large field models 

In this section, we describe the results obtained for large field models. The total number 
of inflationary parameters is now three: the power of the potential p, the potential 
scale M, which is uniquely determined from P*, and the reheating parameter R. This 
accounts for an overall number of parameters of seven given the four base cosmological 
parameters. The best fit model has x 2 — 11252.3. In figure IT?) we have represented the 
(one-dimensional) marginalised posterior probability distributions and mean likelihoods 
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Figure 17. Marginalised posterior probability distributions (solid lines) and mean 
likelihood (dotted lines) for our ACDM base cosmological parameters together with the 
cosmological constant and the Hubble parameter obtained from the exact integration 
of large field models power spectra. 

for the base cosmological parameters. As one may compare with the ones derived in 
the slow-roll models (see figure |2J), their determination is robust and their values remain 
standard. This is just the consequence that, for large field models, there exists values of 
the primordial parameters leading to a power spectrum which allows a good fit of the 
data. 

Figure El shows the one-dimensional posterior probabilities for the primordial 
parameters characterising the large field models. The corresponding two-dimensional 
plots are presented in figure HH1 Let us recall that the quantities p, In R and P* are 
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Figure 18. Probability distributions for the inflationary parameters characterising 
the large field models and the subsequent reheating period. The dotted black 
lines represent the mean likelihood, the solid black lines the marginalised posteriors 
associated with a flat prior on the index p, and the blue dashed lines are the 
marginalised posteriors coming from a flat prior on log p. 



directly sampled by the MCMC whereas M, lni? m d, p en d and a /a en d are derived 
parameters. The curves represent the mean likelihood (dotted black line) and the 
marginalised probability (solid black line) obtained under a flat prior choice on p in 
[0.1, 10], as well as an uniform prior on logp in [—1, 1] (dashed line). 

One sees that models with a potential power p > 4 are now strongly incompatible 
with the data. For values of p slightly greater than, say, 3 — 4, this is due a too high 
level of gravitational waves since the corresponding value of the spectral index n s is still 
compatible with the data. For values of p much larger than 3 — 4, the power spectra are 
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Figure 19. Correlations between various pairs of inflationary parameters 
characterising the large field models and the subsequent reheating period. The coloured 
shaded regions represent the mean likelihoods while the solid contours trace the one 
and two-sigma confidence intervals of the two-dimensional marginalised posteriors. 

too red tilted and their associated n s becomes unacceptable. At two-sigma level, one 
gets the upper marginalised bound 

p<3.1. (159) 

The limiting case p = 1/10 of our prior choice, i.e. an almost Harrison-Zeldovitch 
scale- invariant power spectrum, is still a non-excluded model |56| ISTj . Although the 
maxima of the likelihood and of the marginalised probability are located around p ~ 1.5 
and p ~ 0.5 (respectively), the model p — 1/10 lies in the one-sigma confidence interval. 
This can also be seen in figure EH The fact that values of p < 2 are favoured can be 
recovered in figure |H1 where the corresponding slow-roll approximation predicts a line 
segment below the frontier t 2 — 2e x (corresponding to p — 2) in the (61,62) plane, 
i.e. deeper into the la confidence interval we have obtained in section rOl 

Let us now interpret the distribution obtained for the potential scale M. The black 
solid line represents its marginalised probability when a flat prior on p is assumed while 
the blue dashed line is obtained with a flat prior on log p. The fact that the two curves 
are significantly different signs that this parameter is poorly constrained by the data. 
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On the other hand, the general trend is clearly the same. The value of the peak can 
be understood from equation (|8ip. In fact, the degeneracies between \og(nM) and the 
index p observed in figure HTH can be reproduced almost exactly using this equation from 
which one obtains 

log (kM) ~ - 1.57 + - \ogp - J log (p + 4iV*) 

+ %[1.7-log{p + 4N.)-logp], (160) 

o 

where iV* ~ 50. 

The distribution of p en d can be understood in the same manner. Let us repeat 
that a theoretically motivated upper bound for this quantity is only p cn d < mj,. We 
observe a sharp drop of the marginalised probability for values of ln(K 4 p en d) greater than 
—20. As expected and already mentioned, this is nothing but literally the constraint 
on the energy scale of inflation: in the slow-roll picture, H/m Pl < 1.3 x 10~ 5 . With 
regards to the peak located at ln(ft 4 p en d) — — 22, it can be understood if one notices 
that p en d — M 4 (0 end /m P1 ) p and uses ffTTj) with the most probable p ~ 1. At two-sigma 
level, one gets 

ln( K 4 p end ) < -21.3. (161) 

Let us turn to lni? and/or lni? ra( j. As can be seen on the one-dimensional or two- 
dimensional plots, these two quantities are not constrained by the WMAP3 data. In 
particular, we see the strong influence of the choice of the p prior. With a flat prior on 
p (solid black line), lni? peaks around zero while with a flat prior on lnp (dashed blue 
line) the distribution is almost flat and cut at the edges of its prior. Let us also remind 
that the prior on In R is not of top-hat shape but given by equation ()150|) which involves 
Pend- The data are simply not accurate enough to constraint the reheating phase. Notice 
also that the long tails in the M and p en d distributions are the result of their correlation 
with In R through N* as can be seen in 



Our modelling of the reheating phase by means of lni? only can be slightly 
improved in the special case of the large field models. Indeed, the approach presented 
in section 12.1.21 turns out to be rigourous for those models with the equation of state 
parameter io reh = P/p given by the potential power dependency p through (jTJJ). If we 
express the reheating temperature (J2H|) in terms of the outputs of our computations, 
one obtains 



V4 T „ 30V 4 l/A ( 3 p y /4 

9* J reh — 1— Pend 



v^F rcnd Vp + 8, 

x exp { — In R - - In (/€ 4 p en d) + - In ( - 

I p — 4 |_ 4 V 4\j 



(162) 



In the space [lni?, ln(K 4 p en d),p], a given value of the reheating temperature defines a 
surface. In fact, the above equation can be worked out explicitly and reads 

P 



hxR 



3p 
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Figure 20. One and two-sigma confidence intervals of the three-dimensional 
marginalised probability distribution in the space [(lnre 4 p en( j),p, lni?]. The three 
surfaces represent the locations of constant reheating temperature gl^ 4 T re h, the values 
of which are indicated in the figure. They intersect along a line for p = 4 corresponding 
to a reheating period with an ultra-relativistic equation of state. 



+ y —— In ( K 4 pcnd) . (163) 

The corresponding surfaces are represented in figure EH for three values of the reheating 
temperature, or equivalently, for three values of p re h indicated in the figure. The la and 
2a confidence intervals of the three-dimensional marginalised probability have also been 
plotted. For convenience, we have also represented in figure a few non-marginalised 
two-dimensional sections associated with several values of p. As can be seen in the 
three-dimensional figure, the greater p, the more squeezed are the lines of constant p re h- 
At that point, several comments are in order. Firstly, one notices that the value p = 
4 plays a particular role. In this case, equation ()163J) becomes lni? = (1/4) ln(/€ 4 p end ) 
and lni? does no longer depend on T re ^: this is why all the surfaces intersect along 
the corresponding line as can be seen in figure EH This is also why there is only one 
line in figure |^ for the panel p = 4. Secondly, in this last figure, one also remarks 
that the lines are not horizontal. According to ()163j) . the slope is (p — l)/(3p) (except 
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Figure 21. Two-dimensional cuts of figure EH along the planes of constant p. 
The various lines correspond to (1/4) ln(K 4 p reri ) equal to —45, —35, —25, —15 and 
— (1/4) ln(K 4 /9 on( j) (bottom to top). For the case p = 4, all the lines are merged because, 
in this case, u> re h = 1/3 and the reheating cannot be distinguished from the subsequent 
radiation dominated era. 



for p reri = Pend) and for p = 4 one recovers the factor 1/4. Thirdly, the fact that 
the case p = 4 plays a special role is physically justified. The equation of state 
during reheating reduces to a state parameter u> reh = 1/3 and, as a matter of fact, 
the reheating phase cannot be distinguished from the subsequent radiation dominated 
phase. Fourthly, and as mentioned before, the slope of the lines is different if p re h = Pend- 
This is because, in this case, we have an extra contribution coming from the term 
In (g¥ 4 KT rc hj = In (g¥ 4 K A p cn dj /4. This is especially visible in the top left panel of 
figure ED where all the lines but the one with p rc h = p en d are horizontal. 

Let us now discuss the properties of such a reheating phase. We see in the two top 
panels in figure |2~T1 that for p = 1 and p = 2, there are no constraint on the reheating 
temperature since all the lines intersect the la confidence interval. For p = 1, the 
central value is lni? ~ —15 and corresponds to w re h = —1/3, gl^nT^ ~3x 10~ 8 and 
iV r eh — 20. For p = 1, p re h is unconstrained in the allowed range between p nuc and p en d. 
For p = 2 the central value corresponds to p re h = Pend and, therefore, to iV ren ~ 0, 
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while p re h = Pnuc is on the lower edge of the one-sigma confidence interval. Although 
no constraints can be put on T re h, these plots show that large field inflationary models 
are compatible with the data in the sense that the accelerated phase and the reheating 
period, which share the same field potential, can indeed coexist without leading to 
contradictions when confronted to the CMB observations. The cases p = 3 and p = 4 
are slightly different since the lines of constant reheating temperature are outside the 
la contour for p = 3 and outside the 2a for p = 4. However, these models were already 
strongly disfavoured because of their "high" value of p. 

4-3. Small field models 

In this subsection, we study the exact numerical integration of small field models. The 
number of inflationary parameters is now four leading to an overall number of eight 
model parameters [see equation ([84)1] . The best fit model corresponds to x 2 — 11252.2, 
a value very similar to the one obtained in the large field case albeit involving the 
additional parameter /i. 

In figure E21 we have plotted the posteriors obtained for the base cosmological 
parameters. As previously, the dotted black curves are the mean likelihoods, the thick 
solid black curves represent the marginalised probabilities obtained under a flat prior on 
the parameter /i/m P1 in the range [0.1, 10] and, finally, the thin solid red lines refer to 
the marginalised probability stemming from a flat prior on /i/m P1 but now in the range 
[0.1, 100]. The most probable values of these parameters are compatible with what was 
found in the slow-roll models and the exact integration of the large field models. It is 
nevertheless interesting to notice that the determination of these distributions is not 
completely insensitive to the choice on the prior. This is particularly the case for the 
base parameters Vl^h 2 and r, which is not really surprising since both of them have 
significant correlations with the spectral index n s . 

In figure |2H we have represented the posterior distributions obtained for the 
inflationary parameters. The conventions are the same as in figure 12*21 The blue dashed 
line is the marginalised probability with a flat prior on log(/i/m P1 ) and is represented 
only when it leads to different results. On the top left panel, one sees that the parameter 
p is not constrained since the corresponding distribution is basically flat. Let us notice 
that a flat prior has been chosen for p in the range [2.4, 10], the case p — > 2 being very 
difficult to handle numerically due to computational accuracy limitation. This comes 
from the exponential behaviour of the e\ parameter as we discussed in section 13.41 
In fact the case p ~ 2 is interesting because one notices in figure ESI that the value 
of the marginalised probability depends on the prior in this regime. When the flat 
prior on fi/m Pl lies in [0.1, 10], then the values p ~ 2 are slightly disfavoured but this 
tendency completely disappears when the prior is extended to [0.1,100]. This effect 
can be physically interpreted by means of the slow-roll approximation presented in 
section 13.41 In the limit fi/m Pl <C 1, the first slow-roll parameter is exponentially 
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Figure 22. Mean likelihoods (dotted black lines), marginalised posterior probability 
distributions from a flat prior on n\x in [0.1, 10] (thick solid black lines) and marginalised 
posterior probability distributions from a flat prior on k[x in [0.1, 100] (thin solid red 
lines) for the ACDM cosmological base parameters. The cosmological constant and 
the Hubble parameter posteriors are also represented. They have been obtained from 
a direct numerical integration of the small field models power spectra. 



suppressed while the second one is given by equation ()102|) . i.e. 




(164) 



In this limit, there is no means to comply with the observational constraints because 
€2 becomes large which is strongly disfavoured by the data. However, as shown before, 
the above formula is incorrect if /i/m P1 is not small. When the slow-roll parameters 
are correctly evaluated, we have already pointed in section 13.41 that a compatible value 
of e 2 can be obtained provided a large enough value of /x/m P1 is chosen. Therefore, if 
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Figure 23. Marginalised posterior probability distributions associated with a flat 
prior on k/j, in [0.1, 10] (thick solid black lines) and a flat prior on k/j, in [0.1, 100] (thin 
solid red lines) for the primordial parameters characterising small field models and 
the subsequent reheating phase. The blue dashed line represents the marginalised 
probabilities with a flat prior on log(K/i) in [—1,1] and is drawn only when it is 
significantly different from the others. The dotted vertical line in the top right panel 
indicates the relevant range of values for the k/i G [0.1, 10] prior. 
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the prior is too narrow, then such a large value of /j,/m P1 cannot be reached and the 
case p ~ 2 is disfavoured. On the other hand, if the fi/m P1 prior allows for large enough 
value, then the tilt of the power spectra can be made perfectly compatible with the data 
by singling out high fi/m P1 values. The case p ~ 2 becomes perfectly compatible and, 
hence, the posterior distribution on p flattens. Therefore, from a data analysis point of 
view, we conclude that the model with p ~ 2 is not excluded at all. From a theoretical 
point of view, the situation is less clear since one could argue that the small field models 
require [i < m Pl . Note that, for p = 2, it is necessary to have [i/m m > 0.25 in order for 
the energy scale M to be above the TeV (but /i/m P1 > 0.03 only for p = 2.1). Finally, 
let us stress that we always have 4>/fi <C 1 and this is mandatory because the small 
field potential (J83jl should be viewed as the first leading terms of a Taylor expansion. 
However, this does not require /i/m P1 <C 1. This condition appears when one requires 
the vev of the inflaton field to be explicitly smaller than the Planck mass: a theoretical 
prior choice but of some interest since we have just shown that it modifies the p posterior 
distribution close to p — 2. The question of knowing whether very large inflaton vev in 
comparison with the Planck mass makes sense is controversial [8$ ■ 

The top right panel in figure ED represents the posterior associated with the 
parameter fi/m Pl . As can be seen, this parameter is basically unconstrained. The 
vertical dotted line is the upper limit of the posterior derived under the flat prior in 
the range [0.1, 10]. In this case the respective posteriors obtained from both the prior 
choices match. We notice that, if one focuses on the range [0.1,10] only, then there 
is a tendency for a "large" k/z, i.e. close to the upper bound K/i ~ 10. As explained 
before, in the slow-roll language, very small values of K[i are indeed disfavoured since 
they would lead too large values of the parameter e 2 . In the full range [0.1, 100] the K/i 
posterior is essentially flat and there is no significant constraint on the values of K[i. 

Let us now turn to the energy scale M and the energy density p en d at the end 
of inflation. A first remark is that these parameters are poorly constrained and their 
corresponding posterior distributions are strongly dependent on the prior choice. Let 
us try to interpret these curves in more details making use our slow-roll understanding. 
If fi/m Pl <C 1, then, in the slow-roll approximation, the potential scale M reads 

-2(p-l)/(p-2) / „ N 2p/(p-2) 

. (165) 



M\ ^^5p 2 Q 2 Tms _ PS \ AT p(p-2) - '(jji 



m P] J 327T T 2 



87T 

As already discussed in section I3~H this energy scale can be very small. For fi/m Pl > 
1, following the discussion after equation (|99|) . one may use an approximation of 
equation ()107|) by using the limit 0* — > 1. One gets ei ~ 0(1) which gives a very 
rough order of magnitude estimate (see figure ITUj) . Finally 
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that is to say log(/tM) ~ —1.9. This is in good agreement with what is observed in 
figure 1221 If we use the prior k\i in [0.1, 100] (thin red solid curve) then we allows 
for larger value of k/i than with the prior k/i in [0.1, 10] (thick solid black curve) and, 
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Figure 24. One and two-sigma confidence contours of the two-dimensional 
marginalised posterior probability distributions (point density) in the plane 
[log(«M),p] (left panel) and [ln(«; 4 pend), In -R ra d] (right panel). The correlations with 
Kfi are indicated by the color bar. 



consequently, the marginalised probability is shifted towards the large energy scales by 
volume effects in the parameter space. If, on the contrary, we use an uniform prior 
on log(/i/m P1 ) (dashed blue line), then we favour the small values of /i/m P1 and the 
marginalised probability opens up in the direction of small M/m P1 , in accordance with 
the previous considerations. These ones are also confirmed by inspecting the left panel 
of figure EU For a fixed value of p, the larger k/j, is, the larger \og(nM) becomes. 
Moreover, when k\i > 1, log(nM) has a very weak dependence on p as one may guess 
from (jlfifij) . Finally, let us stress that the same interpretation also applies to p en d- In 
addition we recover the same cutoff at ln(K 4 p en d) — —20 as for large field models, which 
is once again simply the manifestation of the constraint on the level of gravitational 
waves, i.e. in the slow-roll language the result of e\ < 0.022. 

Finally, let us discuss the quantities related to the reheating phase, in particular 
In R. To start with, it is important to remark that contrary to the large field models, the 
quantities characterising the reheating phase cannot be related to the free parameters 
of the small field potential. In particular, the formula (j!4j) is specific to monomial 
potentials and can no longer be applied to the present case. Now, from figure ESI we see 
that the quantity InR has a lower bound slightly above the lower prior bound and seems 
to be constrained. On the contrary, the upper limit only comes from the upper bound 
of the prior. Let us recall that the lower limit of the prior comes from the requirement 
Preh > Pnuc which approximately means lni? > —46. Here, we obtain the marginalised 
lower bound 



lni?>-34, (167) 
at two-sigma level of confidence (and InR > —22 at one-sigma). These values are 



Inflation after WMAP3 



58 




Figure 25. One and two-sigma confidence intervals (solid contours) of the two- 
dimensional marginalised posteriors (point density) in the plane [lni?, ln(/c p e nd)]- The 
colour code is the same as in figure El The reheating predictions are represented in 
the particular cases where the state parameter u> ro h = — 1/3 (top left panel), u> ro h = 
(top right panel), u> re h = 1/3 (bottom left panel) and iu rc h = 1 (bottom right panel). 
The various solid lines correspond to different values of the reheating temperature 
(l/4)ln(K 4 p roh ) ranging from -45, -35, -25, -15 to -(1/4) ln^Vnd)- 



obtained under the prior k/i in [0.1,100] which weakens the lower bound of lni? 
(see figure |23). The corresponding bounds for the derived parameter i? rad reads 
lni? rac j > —31 at 2a. Note that these constraints are recovered on the shape of the 
ln(a /a e nd) posterior at large values. It is important to notice that these limits are 
modified when the prior is changed, reducing the upper bound on the Kfi prior would 
induce tighter constraints on lni?. The small field models are the only class of models 
for which we find some non trivial, albeit weak, constraints on the reheating parameter. 
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Table 1. Theoretical predicted slow-roll values of lni? for small field models with 
p = 4 as a function of £2 and for the mode k/a^ = 0.05 Mpc -1 . Some of these values 
are not compatible with the data (see text). 



Let us now try to understand the physical origin of this bound in more details. At 
first glance it may be surprising to get a constraint on a parameter which is not directly 
involved in the generation of the primordial cosmological perturbations. In fact, as we 
show below, this is a side-effect coming from the constraint existing on the scalar power 
spectrum tilt, that is to say the on the spectral index n s — 1, or e 2 in the small field 
slow-roll approximation. 

In the case of small field models, the first slow-roll parameter is extremely small. 
Using equation ()146|) . one gets 



lni? ~ -58 + In 



A (Mpc" 1 ) 
a 



N k -\\n^, (168) 

^ ^end 



where k/ao is a given physical scale, N k the number of e-folds between the time at 
which the scale k leaves the Hubble radius during inflation and the end of inflation, 
and V k denotes the value of the inflaton potential at Hubble exit. This equation means 
that lni? is directly related to the time of Hubble exit. For different values of lni?, 
the modes of cosmological interest today probe different ranges of field value along the 
inflationary potential. On the other hand, not all the field values are compatible with 
the data, simply because in a region where the potential is too steep (or even too flat), 
the spectral index may take unacceptable values when compared to the observations. 

One can quantitatively test the previous discussing. In practice, we choose 
k/a = 0.05 Mpc" 1 , and one may approximate V k by M 4 . Under these assumptions, 
the effect described before only comes from the values taken by N k . The slow- roll value 
of 4>end is known from ei(0 e nd) — 1 whereas the value of <p k is obtained by solving the 
equation e 2 (</>&) = e 2 | b s - Recall that the function e 2 (</>) is explicitly given in (}85|) and 
£2 |obs is supposed to be a fiducial value of the second slow-roll parameter compatible 
with the data. Using (J9"3~j) and (}9"4"j) for end and (fi k uniquely determines N k and thus 
lni? by means of ()168j) . Since only a limited range of e2| bs values is compatible with 
the data, one indeed expect some constraints on the reheating parameter lni?. For 
instance, all lni? values associated with e 2 | b s = 0.07 should be on the edge of the lni? 
marginalised posterior. 

In table [TJ we have derived the expected values of lni? with respect to e 2 | bs in the 
pure slow-roll approximation for p = 4 and for different values of Kfi. For nfi = 15, 
the value e 2 = 0.07 is on the edge of the confidence interval derived in section 13.21 
and corresponds to lni? ~ —40: this is compatible with the drop of the lni? posterior 
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Figure 26. Marginalised posterior probability distribution for the reheating energy 
scale ln(ftPre;h) assuming a constant state parameter w r ch- For the extreme case 
Wrch =— — 1/3, the small field models prefer a reheating temperature T re h > 2TeV at 
95% confidence level. 



distribution around the same value seen in figure I2B1 It is interesting to note that for 
smaller values of k/i, the relevant iV fc and lni? are pushed toward higher values (see 
table GJ. This confirms that, if one has the theoretical prejudice that fi/m Pl < 1, then 
the bound on In R is much tighter. Finally, to conclude this discussion, one can estimate 
how much the change of the Kfi prior propagates to the other derived parameters. From 
equation ()146|) 

A (ln^-) =-A(]ni2) + jA[ln(« 4 Pend )], (169) 

and from table d and figure ESI we see that when one changes the prior from k/j, in 
[0.1, 10] to Kfi in [0.1, 100], one has A (InR) ~ -10 and A [In (K 4 pend)] ~ +4. Therefore, 
one gets A [ln(a /a e nd)] — +12 which is also the variation observed in the bottom left 
panel in figure 1231 

Let us now try to propagate the weak constraint on In R to the phase of reheating 
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itself. The reheating temperature can be evaluated as in (J162j) except that, as already 
mentioned, there is no longer a link between the state parameter w re ^ and p. However, 
for a constant u> re h, one still has 

1 I 1/4 m \ 3 + 3w re h , D 1 + 3u> re h , / 4 \ 

In [gj KT veh ) = - InR- - — Ib(k p en d) 

1 - 3w rch 2 - 6w Teh 

I 1 ln ( 3 + 3uJ ™A | ln 3Ql/4 (170 ) 
1 - 3u> re h \5-3w veh J 

For a fixed value of T re h (and of u> re h), equation ()17()|) represents a line in the plane 

[lni?, ln(K 4 p end )]. These lines are plotted in figure EH for four values of u> re h and 

for the reheating temperature (1/4) ln(/t 4 p reh ) equal to —45, —35, —25, —15, and 

— (1/4) ln(«; 4 pend)- As can be seen, the bound InR > —34 does not tighten significantly 

the allowed T re h values. The only lower limit that can be extracted concerns the 

rather extreme reheating limit w re h — —1/3 for which the lines corresponding to 

gy 4 T reh < 1 MeV lie outside the two-sigma level contour over \nR and p en d. The 

marginalised posteriors on T re h can be derived by importance sampling ^l] from the 

previous analysis and are plotted in figure EEl for four values of the parameter w re h. 

Notice that we have avoided the pure radiation-like case w re h = 1/3 which is problematic 

in the inversion formula ()17()|) . as previously discussed. In the extreme case w re h — —1/3, 

the reheating temperature satisfies 

^/ 4 T reh >2TeV, (171) 

at 95% of confidence. Of course, such a bound close to the TeV scale is already disfavored 
by BBN and may therefore seem irrelevant. However, we would like to point out that, 
in the present context, it was obtained from the CMB data only and to our knowledge, 
this is the first time that this can be done. 

Let us conclude by noticing that, with better data (for instance coming from the 
Planck satellite jEH]), it is likely to obtain much better constraint on the reheating 
temperature in the framework of small field models. As shown above, this is the result 
of the influence of In R on the observed value of the the spectral index. Therefore, much 
tighter constraints on n s would probably lead to a narrower range of allowed reheating 
temperatures. 



4-4- Hybrid Models 

Let us now turn to the exact numerical integration of hybrid models. This class of 
models involves nine parameters in total, the four base cosmological parameters plus 
five inflationary parameters. The best fit model has a x 2 — 11257.2, a value which 
is larger than the one obtained for large and small field models, the difference being 
Ax 2 — +5 in spite of one extra parameter with respect to small field models and two 
with respect to large field models [see equation (jlU9j) ]. It is therefore fair to say that 
this class of models does not fit the data as well as the two previous ones. 

In Fig. |2Zl we have represented the posterior distributions for the base cosmological 
parameters. The most striking features of these plots are the values of the baryon 
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Figure 27. Marginalised posterior probability distributions (solid black lines) and 
mean likelihoods (dotted lines) for the ACDM base cosmological parameters obtained 
for the hybrid models. The cosmological constant and the Hubble parameter derived 
posteriors are also shown. 

number density Q^h 2 ~ 0.024, and of the optical depth, r ~ 0.14. They should 
be compared to what was obtained before, say for large or small fields model where 
Q^h 2 ~ 0.0225 and r ~ 0.1. In particular, such a large value of Q\> exacerbates the 
tension between the results obtained with the CMB and the BBN. Indeed, one has [HHl 
0.020 < Vl h h 2 < 0.024 from the deuterium measurements and 0.007 < Vt h h 2 < 0.014 
from 7 Li. Even if one considers only the Deuterium measurements, then the previous 
result pushes fl^h 2 towards the upper limit of the BBN allowed range. 

As discussed in section ET51 in the slow-roll context, the hybrid models are associated 
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with a blue tilted scalar power spectrum, i.e. typically n s > 1. Although not totally 
excluded by the WMAP data, this situation is nevertheless disfavoured. Since the 
parameter n s is degenerated with the two parameters fl^h 2 and r, the more likely way 
to accommodate such primordial power spectra with the observed CMB anisotropics is 
to compensate the n s > 1 values of the spectral index by modifying these two quantities. 
In particular, in order to compensate an excess of power at small scales (blue spectrum), 
one has to increase Q^h 2 and r as observed. 

In figure "2*""". we have plotted the posteriors obtained for the inflationary parameters. 
As already discussed, the hybrid models involve an extra parameter, namely the value 
of the inflaton field at which the instability occurs and inflation stops. Also, we have 
to pay attention to probe only the vacuum dominated regime which is, by definition, 
what we mean by hybrid inflation. Otherwise, the situation would be very similar to 
the large field models studied before. 

From the slow-roll analysis of section 13*31 when /i/m P1 < p/(8i/7r), the field has 
to satisfy < 07 (see also figure [131). However, for a chosen value of cri too close to 
4>~ , it is not guaranteed to obtains a total number of e- folds larger than the required 
60. In the following, we define 0i 2 o to be the maximum value of cr i leading to at least 
120 e-folds of inflation (0 is decreasing as hybrid inflation proceeds). In other words, 
cri should be smaller than 12O . 

If /i/m P1 is large enough, then <p~ does no longer exist. In this case, 0i 2 o is defined 
from the value of corresponding to the inflexion point of the potential, i.e. to e 2 = in 
the slow-roll framework, which is also the field value maximising t\. This is a convenient 
way to identify the vacuum dominated regime for several reasons. Firstly, e 2 = 
corresponds to 0//x = 1. Therefore, the regime e 2 < means that 0//i < 1 and, hence, 
the term M 4 dominates in the hybrid potential. Secondly, e 2 < means that the kinetic 
energy is decreasing absolutely, and relatively to the total energy which, again, is exactly 
what we expect in a vacuum dominated regime. In the following, cr ; will be always 
measured in units of 0i 2O . 

As seen in figure EHl the values of the hybrid inflationary parameters are not really 
constrained. The posterior distributions for p or for the k/i are basically flat although 
small values of Kfi are slightly disfavoured. In the same manner, the critical value cr i 
is not constrained but values close to cr j ~ 0i 2 o seem to be slightly favoured. With 
regards to the energy scale M and the energy density at the end of inflation, p cnd , 
we observe the usual cutoff for the large values of these quantities, as it was the case 
for large and small field models. This is still due to the constraint on the amount 
of gravitational waves and in terms of slow-roll, on e±. These quantities are poorly 
determined as they strongly depend on the choice of the prior on log0 cri . In these plots, 
the solid black lines represent the marginalised posterior probabilities obtained with a 
flat prior on log(0 cr i/0i 2 o), whereas the dashed blue lines correspond to a flat prior on 
cr i. The order of magnitude of the results can be very simply estimated. For instance, 
for the scale M, since we are in the vacuum dominated regime, the term (0//x) p can be 
neglected and V ~ M 4 . Then, the WMAP normalisation leads to the same equation as 



Inflation after WMAP3 



64 





Figure 28. Marginalised posterior probability distributions (solid black lines) and 
mean likelihoods (dotted curves) of the hybrid model inflationary parameters drawn 
from a flat prior on log<^ cr j in the range [—2,0]. The blue dashed lines represents the 
marginalised posteriors obtained from an uniform prior on </> cr j whenever they do not 
match with the former ones. 
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Figure 29. Two-dimensional marginalised probability (point density) in the plane 
[log(fcM),f>] for various values of cr j measured in units of 4>i2a and indicated by the 
colour bar. The one-sigma and two-sigma confidence interval are represented by the 
solid contours. 

for small field models, namely (jlfifij) . which implies M/m Pl ~ 10~ 4 — 10~ 2 , in agreement 
with figure EEl Concerning the reheating parameter, there is no constraint at all. The 
distribution of In R is flat and is only cut by the prior. This is of course the same for the 
posterior distribution of the derived parameter ln(a /a en d). Therefore, for the hybrid 
models, nothing can be said on the reheating with the WMAP3 data. 

Finally, figure |2Hl represents the two-dimensional marginalised probability in the 
plane [log(ftM),p] for various values of <fi cri measured in units of 0i 2 o- We see that the 
allowed range for M is tighter for smaller values of p and that models with large values 
°f ^cri/0120 tend to have a scale M larger than models with smaller cr i/0i2o- Such 
an effect is expected since the potential scale M , in the vacuum dominated regime, is 
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Figure 30. Marginalised posteriors (solid curves) and mean likelihoods (dotted) for 
the cosmological parameters in the running mass inflation model. 



roughly given by 

m P1 / 2 1 z 

Now, if we decrease cr i then we decrease 0* and hence the observed value of e\ according 
to figure Uni (bottom left panel). As a result, decreasing cri implies decreasing the energy 
scale M, as observed in figure EH1 
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4-5. Running mass models 

Finally, we now turn to the numerical integration of the running mass models for which 
the best fit parameters lead to x 2 — 11252.3. As discussed in section I3~T)| these models 
are characterised by three additional parameters, c, 0o an d CT i m comparison with the 
large field models. The total number of inflationary parameters is therefore five, as 
for the hybrid models, accounting for an overall number of nine parameters. Recall 
however that cr j comes from the choice to stop inflation by instability and according 
to the different versions of this inflationary scenario, this is not always necessary. 

Our numerical integration of these models appears to be limited in the parameter 
space due to finite accuracy issues. As pointed in section l3~H| the slow-roll approximation 
gives a field evolution 4>(N) involving a double exponential behaviour with respect to the 
number of e-folds [see equation ()132|) ]. As a result, integrating the equations of motion 
along 60 e-folds may require a computing accuracy much smaller than the quadruple 
precision computing bound 10 -32 . In order to satisfy this computing requirement, we 
have considered a rather limited range for the MCMC running mass parameters: an 
uniform prior has been chosen for log(ft0o) in [—3,0], for c/2 in [—0.01,0.1] and for 
log(/t0 cri ) in [-1.8,0]. 

It may be convenient to compare this approach to the one used in the literature jH21 • 
In that reference, when concerns with the leading order only, the slow-roll equations 
n s — 1 ~ 2(s — c) and a s ~ 2sc are assumed. Motivated by theoretical priors, an 
uniform sampling on s and c in [—0.2, 0.2] is performed and used to compare the model 
predictions with the data, presumably by assuming a power-law power spectrum whose 
tilt and running are uniquely determined from the s and c values. An advantage of this 
method is that one can consider a wide range of values for the parameters c and s. Our 
method differs from the fact that no assumption is made on the shape of the primordial 
power spectrum. Therefore, the exact numerical integration requires a sampling on the 
fundamental model parameters c, O > 0cri- The price to pay, as mentioned before, is 
that the resulting domains probed by the derived parameters s and c are reduced to 
maintain the required computational accuracy. 

The one-dimensional marginalised posteriors for the base cosmological parameters 
are plotted in figure EH while the inflationary parameters are represented in figure 
As for the other inflationary models which reasonably fit well the data, the cosmological 
parameters are centred to their fiducial values. Concerning the primordial parameters, 
we recover the standard bounds on the energy scale of inflation and the power spectrum 
amplitude, while the reheating parameters are not constrained. On the other hand, the 
running mass parameters exhibit distorted distributions and large values of 0o seem to 
be favoured. A priori, the posterior on c/2 could be associated with some confidence 
intervals. However, as we show in the following, the patterns observed in figure I3T1 for 
the running mass parameters are mainly dominated by the prior choices and the multiple 
degeneracies that arise between the model parameters and their observable effects. This 
means that these parameters cannot be robustly constrained in a prior independent way 
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Figure 31. Marginalised posteriors (solid curves) and mean likelihoods (dotted) for 
the running mass inflationary parameters. Apart from the power spectrum amplitude 
P* and its derived parameters, these posteriors are poorly constrained by the data and 
mainly result from the prior choice (see text). 



with the current data, at least in our prior range. 

The degeneracies between the parameters are clearly drawn in the two-dimensional 
probability distribution (point density) plotted in figure E21 The correlations between 
c/2, log(0 C ri) and log(0o) appear as the preferred "islands" and "directions" explored 
by the MCMC. The topology of these correlations can be understood from the slow-roll 
approximation detailed in section ETTA Indeed, since e\ is generically small for this class 
of models, equation (J138J) implies 1 — n s ~ e 2 . Since the spectral index n s is constrained 




Figure 32. Two-dimensional marginalised probability (dot density) in the plane 
[log(l/</> ), c/2] for various values of </> cr ; (colour bar). The one and two-sigma 
confidence intervals are the solid contours. 



by the data, one can recover the range of model parameters that lead to such an observed 
value of n s , or €2 in the present case. 

From the expression of s derived in (jl37J) . a given value of €2 defines a surface in 
the three-dimensional volume [log(/t0 cr i), c/2, \og(l / n<fio)} given by 



log(K0 cri ) = e 



-cN. 



(173) 



4 \2J V^o, 

for a given AT*. This equation explicitly gives the degeneracies between the running mass 
parameters leading to the same €2 values. However, it is only an approximate solution of 
the slow-roll equations valid as long as k0 1 and some precautions should be in order 
when dealing with the RM4 and RM2 models. In fact, one may alternatively invert ()130j) 
numerically for a given iV* and use the equation (j!28|) to recover the surfaces of constant 
€2 ■ These surfaces are represented in figure ESI for three values of €2 corresponding to 
the mean value (e2 = 0.034) and to the two-sigma confidence bounds (e^ = —0.029 and 
t2 = 0.074) obtained for the second order slow-roll expansion and under the Jeffreys' 
prior on e\. If we compare figure E3 with figure ESI we recover that the MCMC spread 
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Figure 33. Isosurfaces of constant €2 in the three-dimensional parameter space of 
the running mass inflation model. The blue surface is associated with &2 = —0.029 
and the dark red one with €2 — 0.074. These two values correspond to the two-sigma 
confidence intervals at next-to-leading order in slow-roll expansion and for the Jeffreys' 
prior on e\. The yellow surface corresponds to £2 = 0.034. The black mesh corresponds 
to <fi — <fio and marks the border between RM1, RM3 on one hand and RM2, RM4 on 
the other hand. Finally, the region c > corresponds to RM1 and RM2 while c < 
contains RM3 and RM4. The value = 40 has been used. 



out along the surfaces associated with e 2 ~ 0.034. 

Firstly, let us notice that the folded yellow surface in the region — log(/«0o) — 4, 
c < corresponds to models with small values of K0 cri and kc/) but (and this cannot be 
seen in the figure) with large values of Kip- in (this is the RM4 model). These models are 
compatible with the data but they may suffer a theoretical inconsistency since they are 
associated with vev larger than the Planck mass as inflation proceeds. Note that these 
parts of the yellow surface cannot be found by using (jl73|) since its domain of validity 
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breaks down in that particular case. 

Figure ESI demonstrates that the constraints appearing on the one and two- 
dimensional marginalised probabilities on the running mass parameters are mainly due 
to a preferential volume selection in the parameter space by the prior choice. Indeed, 
if one would extend the range of log(K0 cri ) below —3, then the dark red surface would 
penetrate the region of higher c values and would render these c values compatible with 
the data. The conclusion is similar for the blue inverted half-pipe around the value 
c = 0. Extending the prior of K0 cr i and k0o would render the region around c = 
compatible with the data. However, the strong variations of K0 cri and k4>q required to 
get a reasonable value of €2 for an arbitrarily chosen c suggest that some amount of 
fine-tuning is required between these parameters. 

Nevertheless, for the running mass models, the main conclusion is that the 
dependence on the priors prevents us to establish reliable constraints on the free 
parameters. To avoid this difficulty, it would be necessary to integrate numerically 
the spectrum for very small values of ncj) cr \ and k0o but, as noticed at the beginning of 
the section, this is a non-trivial technical issue and some theoretical lower bound should 
also be set. On the other hand, one can say that these models remain compatible with 
WMAP3 data. 

5. Are the power spectra really featureless? 

5.1. Basic equations 

Recently, the possibility that the power spectra could contain non-expected features 
has been widely discussed. This question arose because of the presence of the so-called 
cosmic variance outliers in the first year WMAP data. Since these litigious points have 
disappeared (at least some of them) in the new data, it seems at first sight that this 
is no longer an interesting issue. This conclusion is also consistent with the analysis of 
the WMAP team [2] which quotes a x 2 drop of A% 2 ~ —4.5 or A% 2 ~ —9.5 according 
to the type of features considered. Notice that, although Ax 2 — —4.5 seems indeed 
not of much interest, Ax 2 — —9.5 deserves at least some attention. In this article, 
we would like to consider this question further and derive the marginalised bound 
the parameters associated with these features satisfy. A tool that has been used in 
order to address this issue is to confront the trans-Planckian power spectra to the 
CMB data. The trans-Planckian problem originates from the fact that, due to the 
exponential expansion during inflation, the scales of astrophysical interest today were 
below the Planck length at the beginning of inflation EBE21I23II21I23I2SIE3I2BII22I 

rmni rnm rrn^i rrnni rrrpi mrKi rmm mT7i rmHi rrrn^i rrmi rrrn rrrsi rrn^i mm rrr^i mm rrrrj . 

In this regime, the framework utilised to derive the inflationary predictions (namely 
quantum field theory in curved space-time) breaks down. This is similar to what 
happens in the context of black hole physics and the derivation of the Hawking 
radiation |1181 11191 H^Uj . If we denote M c the physical scale at which new physical 
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effects are supposed to become relevant, one finds that superimposed oscillations appear 
in the power spectra the amplitude of which is controlled by the H to M c ratio, the 
only two scales available in the theory ||. These superimposed oscillations are then tested 
against the data and used as a tool to detect non-trivial features in the primordial power 

spectrum H2SJ EE22 U2SI U22I U3DI E3U EE32 ESI rTMfTTTKj . 

It was shown in Ref. |107| 11291 H31| H3Uj that a crucial point is the fact that the 
amplitude and the frequency of the oscillations are independent quantities. Since this 
question was sometimes not fairly appreciated in the recent literature, we would like to 
briefly review where this comes from. In particular, we argue that considering dependent 
amplitude and frequency is theoretically not justified in the framework of the so-called 
"minimal approach" . Moreover, postulating that the amplitude is inversely proportional 
to the frequency is phenomenologically restrictive since this has the consequence that 
the region of interest in the parameter space is missed. 

In the minimal approach, only the initial conditions for the primordial perturbations 
are modified while their equation of motion is left unchanged (this would no longer 
be the case if, for instance, we had modelled the new physical effects by a modified 
dispersion relation). The initial conditions are not the standard ones because they are 
fixed when the wavelength of a given Fourier mode becomes equal to a new fundamental 
characteristic scale l c . The time r\ k of mode "appearance" associated with a comoving 

wavenumber k stems from the condition 

2tt 2tt 
H^) = T a( Vk )=i c ^—, (174) 

which implies that rj k is a function of k. Therefore, this is different from the standard 
inflationary calculations where the initial time is taken to be r] k = —00 for any 
Fourier mode k and where, in a certain sense, the initial time does not depend on 
k. This additional ^-dependence is at the origin of the appearance of the superimposed 
oscillations. Another crucial question is in which state the Fourier mode is placed at 
the time 77^. At this point, one would like to be as general as possible and we take 

VsAVk) = =F 7= . , , (175) 

^2u StT (r]k) m Pi 

, 1 x , ■ I UsAVk) 4y/n(ck - d k ) . . 

/VW ) =± «V o • 176 

m P1 



The coefficients c k and d k are a priori two arbitrary complex numbers satisfying the 
condition |cfc| 2 — |c4| 2 = 1- Without restricting the physical content of the problem, and 
given the fact that, in the limit M c — > +00, one must recover the standard limit, i.e. the 
Bunch-Davies vacuum, one typically expects that 

c fc = 1 + y k o- + ... , d k = x k o + ... , (177) 

where <To = H/M c . Any other particular choice would be a strong assumption to be 
justified. The only simplification that one may consider is to assume that x k ~ x 

|| Let us notice that superimposed oscillations could also originate from other physical mechanisms [1211 

inainsiinainsi 
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and yk ~ y, that is to say that these two coefficients are not strongly scale dependent 
in the range of scales under consideration. In the following, the parameters x and y 
are considered as free parameters that are not fixed by any existing well-established 
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1. This implies that 



theories except that, of course, they should satisfy \cy.\' 
the amplitude and the frequencies of the superimposed oscillations are independent 
quantities as announced. Finally, the form of the power spectra stemming from the 
previous considerations read, for the scalar modes |108t 1129] 
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and for the gravitational waves 
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(179) 

We see that the 



where ip is the argument of the complex number x, i.e. x 
new power spectra depend, at most, on three new independent parameters, namely the 
amplitude, the frequency and the phase of the superimposed oscillations. In particular, 
their wavelength can be expressed as 

^ = W (180 ) 
k ei 

The derivation of the trans-Planckian corrections in the power spectra ()178|) and 
()179j) assumes that the back-reaction effects are not too important. For consistency, 
the energy density of the perturbations must be smaller or equal than that of the 

\x\ 



inflationary background. This leads to the condition \x\ < v37rm P1 /M c , an estimate 
which is in agreement with the one derived in |l()7j . In order to put numbers on the 
above constraint, we can use equations (jp^j) and (JHZj) together with M c = H/a to arrive 
at 



\x\go < 10 x 



<7n 



;isi) 



where, in order to derive an order of magnitude estimate, the unimportant factors of 
order one have been neglected. It is important to emphasise that the above constraint 
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is only a sufficient condition, but by no means, unless proven otherwise, a necessary 
condition for the validity of the power spectra calculations |llUj . 

5.2. WMAP constraints on the oscillatory parameters 

In order to test the viability of superimposed oscillations in the primordial power spectra, 
we perform an exploration of the primordial parameter space by using MCMC methods 
implemented in COSMOMC , given the third year WMAP data. Our analysis proceeds in 
two steps. In a first part, the trans-Planckian power spectra in f!178p and ()179|) are used 
to seed the CMB anisotropies in the framework of the first order slow-roll expansion 
(see section 13.2. 1J) . In a second part, we reiterate the analysis by using a power law 
primordial power spectrum supporting a wider class of superimposed oscillations, namely 
oscillating with a k/k* power-law dependence. As detailed in 1 1 201 11301 1131) . we use 
a modified version of CAMB to compute the CMB anisotropies whose required accuracy 
and computational time renders impossible an exploration of the full parameter space. 
Along the lines drawn in those references, only the primordial parameter space is probed 
while the base cosmological parameters remain fixed to their best fit values obtained 
from an MCMC analysis without superimposed oscillations. 

5.2.1. Trans-Planckian power spectra As described in the previous section, in addition 
to the first order Hubble- flow parameters e\ and e2, the scalar power spectrum amplitude 
at the pivot scale P*, considering the trans-Planckian power spectra (j!78|) and (|179|) 
accounts for three more primordial parameters. We have chosen uniform priors on the 
overall oscillatory phase tp in [0, 27r], defined as 

$ = —(l + e x )+ip, (182) 

as well as on the parameter \x\o~q in [0, 0.45]. Moreover, in order to sample directly from 
the oscillation frequency, an uniform prior has been chosen on the parameter log(ei/cr ) 
in [1, 2.6]. The other primordial parameters are sampled according to the uniform prior 
choice on ei described in section 13.2. 11 

The converged posteriors for the primordial parameters have been plotted in 
figure 0U As previously mentioned, they have been obtained for a set of fixed 
cosmological parameters tt h h 2 = 0.021, n dm h 2 = 0.0159, 6 = 1.0393 and r = 0.0942 
(implying Hq ~ 72km/s/Mpc) and the MCMC exploration has been stopped after 
approximately 200000 elements for which the generalised Gelman and Rubin /^-statistics 
implemented in COSMOMC |11[ I13b'j is less than 10%. The constraints obtained on 
the standard primordial parameters e\, 62 and P* are stronger than those derived in 
section 13.2. 1| as expected since the cosmological parameters have been fixed to their 
best fit values. 

The overall constraint on trans-Planckian superimposed oscillations is given by the 
\x\o~q marginalised posterior. The vanishing value of this quantity corresponds to the 
standard first order slow-roll primordial power spectra without oscillation and is still 
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Figure 34. Marginalised posterior probability distributions (solid right stairs) and 
mean likelihood (shaded bars) for the trans-Planckian primordial parameters. Note 
that these posteriors are derived under a fix set of cosmological parameters. 

the favoured model given the third year WMAP data. However, as it was the case 
with the first year data the mean likelihood is peaked over non- vanishing values 

of \x\a showing that superimposed oscillations provide a better fit to the data. As 
detailed in the marginalised posterior remains peaked around vanishing values 

due to volume effects in the parameter space: the best fit region occupies a rather 
limited volume that does not take over the accessible volume associated with the non- 
oscillatory models in spite of their lower likelihoods. Nevertheless, the statistical weight 
coming from the highest likelihood values broadens out the posterior distribution on 
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Figure 35. In the left panel, smoothed one and two-sigma contours (solid line) of 
the two-dimensional marginalised posterior probability distributions (coloured dots) 
in the plane [log(ei/ero), M^o]- The highest probable frequencies clearly appear as the 
vertical coloured dot alignments. The right panel represents the cumulative residual \ 2 
between the trans-Planckian best fit model and the vanilla first order slow-roll model. 



\x\(Jo and at two-sigma level one has 

|x|a < 0.38. (183) 

This number has to be compared with the limit derived in |131| . namely |a:|<To < 0.11. 
The best fit volume confinement in the parameter space may be understood on the 
marginalised posterior associated with the frequency parameter log(ei/cro)- As can be 
seen in figure both the mean likelihood and marginalised probability exhibits narrow 
peaks on particular frequencies only thereby leading to a discrete set of best fit sub- 
manifolds in the parameter space. From the WMAP team likelihood code PQElElll], 
the best fit associated with the highest resonance peak leads to an overall x 2 — 11239.9, 
which is a fit improvement of A% 2 ~ — 12 with respect to the first order vanilla slow-roll 
model, for three additional parameters (see figure l33|) . Note that this value is larger 
that the one reported by the WMAP team [2], possibly due to the localisation of the 
highest likelihood resonance in a rather high frequency region log(ei/o"o) — 170. 

Eventually, compared to the first year WMAP data, we still find no evidence for 
superimposed oscillations in the WMAP third year data. However, considering trans- 
Planckian-like superimposed oscillations still significantly improves the fit to the data, 
and more importantly the overall statistical weight associated with these best fit regions 
has increased. 

Let us now study what the above results imply for trans-Planckian physics. As 
already mentioned, the best fit is obtained for log(ei/cro) — 2.23, e\ ~ 2.1 x 10~ 3 and 
\x\o~o ~ 0.268. This implies a ~ 1.2 x 10~ 5 . Moreover, if one uses the value of ei for the 
best fit (this does not mean, of course, that we have detected a non-vanishing e\ since this 
would imply a detection of primordial gravitational waves), then one can estimate M c 
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which reads M c ~ 0.3m P1 . Some serious problems show up when one tries to see whether 
the best fit suffers from a back-reaction problem. Using equation (|181|h one finds that 
this is not the case provided \x\o~o < 3.3 x 1CT 5 . This limit is thus largely violated 
by the best fit. At this point, several remarks are in order. Firstly, strictly speaking, 
this results clearly invalidates the perturbative framework used in order the derived the 
power spectra with the superimposed oscillations. Secondly, as discussed in jl lUf. 1117] . 
the back-reaction is not necessarily a problem as its effect might just "renormalise" the 
vacuum energy during inflation. In other words, it may not necessarily prevent inflation 
to proceed. Thirdly, the presence of superimposed oscillations is not necessarily linked 
to trans-Planckian effects. In that case, the limit given by equation ()181j) simply does 
not apply. But, then, it becomes more difficult to physically motivate the logarithmic 
shape of the oscillations which turns out to be favoured by the data (see below). 

To conclude this section we would like to describe a few intriguing features, although 
not statistically significant, that appears in the previous analysis. As shown in |137j . 
it is worth stressing that multiple resonances in the likelihood are precisely expected in 
presence of an oscillatory signal due to "frequency beating" between the data and the 
model tested. These degeneracies produce a multi-valued function in the recovery of a 
primordial oscillation frequency, but also open a window on the a priori unobservable 
high frequency signals through their lower frequency resonances. Compared to the first 
year data, we find the appearance of new favoured frequencies that can be seen as the 
peaks in the log(ei/cr ) posterior of figure 13^1 or through the dot alignments in figure 1331 
Another intriguing feature concerns the marginalised posterior and the mean likelihood 
of the phase parameter ip. As can be seen in figure H^J values around ip ~ and 
ip ~ 3 are slightly favoured by the data (once again, not in a statistically significant 
way). As shown in |129j . in presence of high frequency oscillations, ip = it modulo 7r 
are the phase values expected to maximise the oscillation amplitude in the multipole 
moments. As already mentioned, even in presence of an oscillatory signal in the data, 
one cannot conclude that it comes from a primordial origin and it may be the result 
of some foreground contamination. However, one would have to explain how such a 
pattern in the inter-multipole correlations arises. 



5.2.2. Power law power spectrum In this section, we reiterate the previous analysis by 
using a phenomenological primordial power spectrum. Only the scalar modes have been 
considered with a power spectrum of the form j!28j 

+ (184) 

In the limit p — > 0, one recovers logarithmic-like dependence in the oscillation frequency. 
The interest of this kind of power spectrum resides in the comparison of models having 
superimposed oscillations with a different frequency dependence in the wave number k. 
In other words, one can test whether the logarithmic oscillations are special or if any 
oscillatory signal can significantly improve the best fit. For the sake of simplicity, we 
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Figure 36. Marginalised posterior probability distributions (solid right stairs) and 
mean likelihood (shaded bars) for the oscillatory power law parameters. Note that 
this posteriors are derived under a fix set of cosmological parameters. 



have not considered a running spectral index and, as in the previous section, only the 
primordial parameter space has been explored once the cosmological parameters have 
been fixed to their best fit value obtained from the fiducial power law scalar power 
spectrum: Q h h 2 = 0.0223, Q dm = 0.1064, 9 = 1.040 and r = 0.0885. An uniform 
prior has been chosen for log(p) in the range [—5, 0.48] (the upper bound corresponds 
to p ~ 3) as well as for log(o;) in [1, 2.5]. The phase ip is sampled from a flat prior in 
[0, 7r] and n s from the range [0.5, 1.5], also with an uniform prior. 

The posterior marginalised distributions given the third year WMAP data are 
plotted in figure EH1 One still observes some particular frequencies improving the 
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fit to the data while their weight on the marginalised probability remains negligible 
(see the posterior). The associated best fit value corresponds to x 2 — 11242.7, 
hence a A% 2 ~ —10 with respect to a fiducial power law model, with however four 
additional parameters. This result is in agreement with the one found by the WMAP 
team 2j. However, it is important to stress that as soon as one consider non- logarithmic 
oscillations, the effect coming from modifying the pivot scale can not longer being 
viewed as a phase redefinition. In the present approach, we have not considered this 
effect and k* has been kept to its standard value 0.05 Mpc -1 . In this respect, the 
posterior on log(p) leads to a two-sigma level upper bound 

p< 0.68, (185) 

in favour of small p value, i.e. logarithmic-like superimposed oscillations. 

6. Discussion and Conclusion 

In this last section, we would like to briefly recap the results obtained in this article. 
In a first step, we have studied the compatibility of inflation with the WMAP3 CMB 
data using the slow-roll approximation. We have found that, at leading order, the first 
slow-roll parameter satisfies at 95% of confidence 

e x < 0.022. (186) 

These constraints implies an upper bound on the contribution of primordial gravitational 
waves, namely rio < 0.21 at 2a. This also leads to an upper bound on the energy scale 
of inflation H/m P1 < 1.3 x 10 -5 . 

The WMAP3 data also constraint the second slow-roll parameter (at the 2a level) 

-0.02 < e 2 < 0.09, -0.07 < e 2 < 0.07, (187) 

the first result being derived with a Jeffrey's prior on e\ while the second is obtained 
with a uniform prior on e±. We see that positive e 2 , hence red spectral index, are slightly 
favoured although a scale-invariant power spectrum remains compatible with the data. 
At the second order in the slow-roll parameter a tendency for e% > is observed but 
this is not statistically significant. Together with e 2 > 0, this would imply a negative 
scalar running. 

Our second step has been to exactly integrate, using numerical methods, the 
inflationary power spectra for four fiducial models. The four models considered were 
the large field, small field, hybrid and running mass scenarios. For large field models 
V((p) oc P , we have found the 2a upper limit 

p<3.1. (188) 

With regards to the subsequent reheating period, constraints on the reheating 
temperature can be found but only for the models that have an index p such that 
they are already excluded. These constraints are therefore not of so much interest and 
are maybe just another indication that these models are now incompatible with the 
CMB data. 
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For small field models, the situation is slightly more complicated. A priori no 
constraint on p and/or /x can be found given the WMAP3 data. But this statement 
is not completely prior independent. Indeed, if one assumes that /x/m P1 < 100, then 
the marginalised probability over \i is flat. But if one considers that fi/m Pl < 10, then 
the case p = 2 is disfavoured. The situation is even more complicated because, if small 
values of fi/m P1 may appear to be appealing from a theoretical point of view, the scale 
H should not be too small in order for the inflation energy scale M to be larger than, 
say, the MeV. Another interesting issue related to the small field models is the reheating 
phase. These models are the only ones for which it is possible to say something on the 
reheat temperature given the WMAP3 data. The bounds are relatively weak since we 
find after marginalisation 



at 95% of confidence. Moreover, they are valid only for quite an extreme equation of 
state during reheating, namely w re h — —1/3. 

We have also studied the hybrid and running mass models. For those models, 
it is even more difficult to say something. Basically, hybrid models are disfavoured 
(Ax 2 = +5) because of their blue spectrum while no prior independent constraint 
can be put on the running mass models because of the strong degeneracy among the 
parameters. But this class of models remain compatible with the CMB data. 

The last section of the paper was devoted to the possible presence of superimposed 
oscillations in the power spectra. The marginalised probability is still centred on 
\x\o~q = 0, i.e. still compatible with no oscillations. At the 2a level, one has obtained 



However, the likelihood is peaked at a non- vanishing value of \x\o~q corresponding to 
Ax 2 = — 12 for 3 extra parameters. This apparent discrepancy is explained by the 
fact that the best fit occupies a small volume in the parameters space. One can 
nevertheless made the following two remarks. First, the overall statistical weight of 
the superimposed oscillations has increased in comparison with WMAP1 data and 
despite the disappearance of some of the cosmic variance outliers. One could have 
expected exactly the opposite. Secondly, we have also tested another functional shape 
for the superimposed oscillations and have shown that it does not improve the fit in 
the same manner, Ax 2 = —10 for 4 extra parameters. This suggests that the precise 
shape of the superimposed oscillations is relevant and that logarithmic oscillations are 
favoured. Notice, however, that the best fit solution, if interpreted in the trans-Planckian 
framework, suffers from a severe backreaction problem. 

Let us conclude this article by a few words about the future. The flow of 
high accuracy data has not yet dried up and the forthcoming CMB experiments, as 
Planck [HH] , will provide us with even more accurate data which will help us to improve 
our constraints on the various inflationary scenarios. Using the slow-roll language, an 
exciting prospect would be to close the contour of the e\ parameter. This would imply 
a detection of primordial gravitational waves and would open the possibility to test the 



T reh > 2TeV, 



(189) 



x\a < 0.38 . 



(190) 
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consistency check r 10 ~ — 5n T , a smoking gun for slow- roll inflation. But even if this 
cannot be done in a close future, the example of small field models has taught us that 
this could also improve our knowledge of the reheating period. Indeed, we have seen 
that the constraint on lni? was directly linked to the constraint on the shape of the 
primordial power spectra. Shrinking the error bars on the tilt could therefore help us to 
put relevant limits on T re h, at least for the small field models but, maybe, also for the 
other class of inflationary scenarios. On the exact integration side, the method we have 
presented could be applied to more complicated models of inflation, especially to those 
transiently violating the slow-roll conditions or to the ones involving several interacting 
fields. Assuming only linear perturbation theory, such an approach directly leads to 
marginalised constraints on the fundamental parameters, as the ones involved in the 
inflaton potential and reheating. A natural extension would be the determination of the 
Bayesian evidence associated with each of the model tested [138, 56J. This would allow 
a statistical meaningful measure to prefer one inflationary model over the others. 
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